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Abstract. An alternative description of quantum scattering processes rests on inhomogeneous terms amended to the 
Schrodinger equation. We detail the structure of sources that give rise to multipole scattering waves of definite angular 
momentum, and introduce pointlike multipole sources as their limiting case. Partial wave theory is recovered for freely 
propagating particles. We obtain novel results for ballistic scattering in an external uniform force field, where we provide 
analytical solutions for both the scattering waves and the integrated particle flux. Our theory directly applies to p-wave 
photodetachment in an electric field. Furthermore, illustrating the effects of extended sources, we predict some properties 
of vortex-bearing atom laser beams outcoupled from a rotating Bose-Einstein condensate under the influence of gravity. 

PACS numbers. 3 . 6 5 . Nk — Scattering theory 
3 . 7 5 . -b — Matter waves 

3 . 7 5 . F i — Phase coherent atomic ensembles; quantum condensation phenomena 
32 . 80 . Gc — Photodetachment of atomic negative ions 



1 Introduction 

The customary approach to elastic quantum scattering phenomena employs a superposition of an incom- 
ing plane wave and an outbound scattering wave that emerges from a localized scattering potential. It is 
sometimes advisable to reduce the complexity of this process by dividing the scattering event into subse- 
quent "absorption" and "emission" stages. The evolution of the emerging wave is then considered separately. 
Obviously, in this description a "reservoir" of particles in the interaction region is required that continu- 
ously feeds the stationary scattering wave. Since the particle number is a conserved quantity in the standard 
quantum picture, we devise a modified approach: In analogy to electrodynamics and other field theories, an 
inhomogeneous "source term" added to the Schrodinger equation allows for particle generation in a finite 
volume. It was Schwinger who introduced the idea of particle sources in field theory. In doing so, Schwinger 
could avoid the use of operator fields. Interestingly enough, he also presented the non-relativistic limit of 



Ref. [BRK97, KBK02] 



such a particle source [Sch73]. Recent examples illustrating the use of the source formalism are presented in 



In our contribution, we inquire into quantum sources that give rise to scattering waves carrying non- 
vanishing angular momentum, which we will denote as multipole waves. In the long wave (or low energy) 
limit, apart from its angular dependence the actual structure of the source becomes insignificant, and the 
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scattering process may be properly modelled using an idealized pointlike source of suitable orbital symme- 
try. Technically, these "multipole sourc es" are o btained from the Dirac (5-distribution (that itself pertains to 
isotropic or s-wave emission, see Ref. [ KBK02 ]) by a simple differentiation procedure outlined in Section ^ 
that grants immediate access to the corresponding multipole wave and currents. 

For scattering waves propagating freely (or in a central potential), multipole sources generate the spher- 
ical waves familiar from partial wave theory (Section Non-trivial res ults emerge, howev er, when the 
scattered particles are subject to acceleration in a homogeneous force field [BBG+98, KBK02]. We present 
analytical expressions for the ensuing ballistic multipole waves and currents that so far have not appeared in 
the literature, and discuss some of their intriguing features in Section |^. These developments are directly ap- 
plicable towards near-threshold photodetachment microscopy, an experimental technique recently introduced 
by Blondel et al. [ BDD96 , BDDV99 1 that allows to observe interference of electron waves on a macroscopic 
scale. Here, we predict the photoelectron distribution in p-wave detachment (Section^. In general, the mul- 
tipole formalism breaks down when the spatial extension of the source becomes comparable to the particle 
wavelength. In the ballistic environment, Gaussian sources provide an important exception since the scat- 
tering waves generated by them apparently converge onto a displaced pointlike "virtual source" |KBK02]. 
Actually, this situation is encountered for an atom laser beam outcoupled from an ideal Bose-Einstein con- 
densate (BEC) that is subsequently accelerated in the earth's gravitational field [MAK+97, BHE99|. Angular 
momentum transfer to the superfluid condensate leads to the formation of vortices that in turn act as sources 
for higher modes of the resulting atom laser The effects of vortices on the beam profile are investigated in 
Section ^. Exact solutions are presented for a single vortex in an otherwise spherically symmetric ideal BEC. 
Furthermore, we calculate the structures imprinted on the atomic beam by the rotating vortex lattices recently 
realized experimentally [ |MCWDO0| , |ASRVK01| , ^CHC02| ]. Finally, some useful mathematical developments 
have been assembled in the appendix. 



2 Multipole sources 

Let us start out with a brief overview of the source formalism and its basic results, emphasizing its kinship 
to conventional scattering theory. Throughout much of this section, we will study the behaviour of idealized 
"multipole sources" that emerge in the long-wave limit. For freely propagating particles, the connection to 
the partial wave formalism is straightforward. 



2.1 A primer on source theory 

In order to moti vate the introduction of quantum sources, we investigate elastic potential scattering as a 
tutorial example [ Mes64 ]. (The source appro ach is re adily extended to more sophisticated problems like the 
atom laser (Section as illustrated in Ref. [ KBK02 ].) In the customary treatment, the total wave function 
■0(r) in the potential V^(r) is decomposed into an incoming wave ■(/'in(r), usually a plane wave, and a scattered 
wave '0sc(r) that may be written as a sum of partial waves '0im(r) of definite spherical symmetry: ijj{r) — 
fpin (j") + tpsc (r) . Obviously, V'in (r ) is not an eigenfunction to the full Hamiltonian H = T + V, but rather to a 
simpler "unperturbed" Hamiltonian TJo = T+U: iJo V^in(r) = Ei/jinir). (Often, one sets U — 0. When long- 



range forces are present, like in Coulomb scattering [ VIes64|, this choice is poor, and U{r) should account 
for the interaction potential. See also Section ^) Consequently, this procedure leads to the introduction 
of the scattering potential W{r) = V{r) — U (r), and the stationary Schrodinger equation reads, as usual, 
[E - H]iP{r) = or: 



[E - Ho] [^i„(r) + Vsc(r)] = W^(r)[V'in(r) + Vsc(r)] 



(1) 
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Since V'in(r) is an eigenfunction to Hq, we may state (jl|) in the equivalent form: 

[E-Ho~ W{r)] ^,e(r) = M^(r)V^i„(r) . 



(2) 



We infer that the scattering wave ijjsc (r) solves a Schrodinger equation for the full Hamiltonian H = Hq + 
W{r), albeit with an additional inhomogeneous term a{r) — W{r)ipin{^) present. 

In view of other inhomogeneous field equations, e. g. Maxwell's equations, the right-hand term (7(r) in 
is identified as a source for the scattering wave V'sc(r). This observation motivates a simple picture for the 
scattering process: The incoming wave ipin{r), via the perturbation W{r), feeds particles into the scattering 
wave V'sc(r) that is governed by the Hamiltonian H. Thus, the decomposition of the wave function into an 
incoming and a scattered part naturally leads to the notion of a quantum source. 

We now turn to the mathematical aspects of (^. Introducing the energy Green function G(r, r'; E) for 
the Hamiltonian H defined via |Eco83]: 



[E-Ho- Wir)] G(r, r'; E) = ,5(r - r') , 



a solution to (Eh in terms of a convolution integral reads: 



^.c(r) = j d\'G(ry-E)a{v') . 



(3) 



(4) 



In general, this result is not unique. However, any two solutions V'sc(r) differ only by an eigenfunction 
'0hom(r) of H: Hip\io„i{r) — £'?/'hom(r). The ambiguity in V'sc(r) is resolved by the demand that G(r, r'; i?) 
presents a retarded solution characterized by outgoing-wave behaviour as r — > oo. Formally, this enforces 
the choice [^S5|]: 



G(r,r';E) = lim 



E-H 



177 



PP 



E -H 



m5{E - H) 



(5) 



where PP(. . . ) denotes the Cauchy principal value of the energy integration. 

Defining the current density in the scattering wave in the usual fashion by j (r ) = ^3 [V'sc (r) * VV'sc (r)] /M 
(where for simplicity we omitted the vector potential A(r), see [KBK01|), the inhomogeneous Schrodinger 
equation (0) gives rise to a modified equation of continuity [BRK97, KBK02]: 



V-j(r) = --3[a(r)*V,e(r)] 



(6) 



Thus, the inhomogeneity (T(r) also acts as a source for the particle current j(r). By integration over the source 
volume, and inserting (^, we obtain a bilinear expression for the total particle current J{E), i. e., the total 
scattering rate: 



J{E) = --3 



J dV J dW(r)*G(r,r';£;)a(r') 



(7) 



Some important identities concerning the total current J{E) are most easily recognized in a formal Dirac 
bra-ket representation. In view of (Bb, we may express J{E) by: 



JiE) = -^5 [{a |G| a)] ^ ^ {a \6iE H)\ a) 



(8) 
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from which the sum rule immediately follows [ KBK02 ] : 



dEJ{E) 



27r 



2n 
li 



dV|cr(r)|^ 



(9) 



(provided this integral exists). In order to connect to the findings of conventional scattering theory, we 
display J{E) in an entirely different, yet wholly equivalent fashion. Employing a complete orthonormal 
set of eigenfunctions \ipfi) of the Hamiltonian H, S{E — H)\ipfi) — 6{E — Eti)\ipfi) follows, and replacing 
|(t) = VFjV'in) we may formally decompose (jsj) into a sum over eigenfunctions; 



(10) 



Thus, Fermi's golden rule is recovered. Another notew orthy con sequence of (^ and (||) emerges in the limit 
of pointlike sources, cr(r) - C6{r - R). We then find [|BRK97|]: 



JiE) = -^\C\'^G{R,-R;E)] = 



(11) 



where n(R; E) = J2fi S{E — Efi)\il!{i(Ti)\'^ is the local density of states of H at the source position R. Equa- 
tion dlT) forms the theoretical basis of the Tersoff-Hamann description of scanning tunneling microscopy 
||TH83i BRK97|]. 



2.2 Multipole sources, waves, and currents 

Evidently, the theory of quantum sources becomes particularly simple for point-like sources. In this case, 
the otherwise bothersome integrations involved in the determination of the scattering wave t/jgc (r) (0) and the 
total current J{E) (Q) become trivial. For the naive choice of a point source, <j{r) — C6{r ~ R), 7/'sc(r) 
is simply proportional to the Green function G(r, R; E) itself, and J{E) follows from the Tersoff-Hamann 
rule ([Tl|). The approximation of a pointlike source is obviously well justified in near-threshold scattering 
{E 0), where the long wavelength of the emerging wave effectively obliterates the internal structure of the 
source. This statement, however, must be taken cum grano salis, for it does not take into account the orbital 
structure of the scattering wave. 

In fact, as discussed below, a simple delta source a{r) ^ S{r — r') invariably leads to a locally isotropic 
emission pattern, i. e., describes scattering into an s-wave. Despite often being appropriate in practice 
[BDD96 |BDDV99 , ECBKOl], conservation of angular momentum may enforce selection rules that restrict 



scattering to higher multipole waves. For these, the idealized point source approach must be suitably mod- 
ified. We proceed by analogy with the well-established multipole formalism in potential theory (or electro- 
statics) as another inhomogeneous field equation. 

Like the scattering wave for a simple point source, the Green function of potential theory equals the 
field created by a source of unit strength located at r', G(r, r') = — l/47r|r — r'|. Additional solutions that 
likewise show a singularity at r' can be constructed via differentiation with respect to the source position 
r' = (x', y' , z'). Of special significance are the multipole potentials $;„j(r, r'): 

, ,, YUer-r') _ K;„,(r-rO 

since they are obviously of (I, m) spherical symmetry. Here, we introduced the harmonic polynomials 
Kim(r) — r^Yijn{f), homogeneous polynomials of order I in the coordinates x,y, z that are eigenfunctions of 
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the angular momentum operator [ ^iob3 1 , MF53 , Mul6€ ]. Interestingly, the sa me polynomial in mome ntum 
space, known as the spherical tensor gradient K/„i(V') = Kim{dx' , dy' , dz') [ Bay78 , Row? 8 , WS83 ], can 
be employed to generate the multipole potentials from the Green function G(r, r'): 



(21 - 1 V 

Kz„(V')G(r,r') = -^—^<^i„,{r,r') 



An 



(13) 



Thus, the spherical tensor gradient imprints the orbital structure onto G(r, r'). Since AG(r, r') = S{r — r') 
holds, we formally obtain from (O): 



A$,™(r,r') = 



47r 



{21 - 1) 



K,,„(V')<5(r-r') 



(14) 



(Note that the differentiation proceeds with respect to r'.) Apart from prefactors, the multipole potentials 
^im{r, r') are thus generated by the spherical tensor gradients of the delta distribution. 

Accordingly, for the purpose of the quantum source problem, we define multipole point sources 6im{r — 
r') via [|Bay78i|Row78t: 



(15) 



Since K;,„(V') and the Hamilton operator H{r, p) always commute, the inhomogeneous Schrodinger equa- 
tion: 



(16) 



is formally solved by the multipole Green functions Gim{r,r' ; E) that are available from G{r,r';E) by 
differentiation: 



GUr,r';E) ^KUV')G{T,r';E) 



(17) 



(We remark in passing that G(r, r'; E) differs from the s-wave multipole Green function Goo(r, r'; E) only 
by a factor VItt.) The multipole point sources (5/„i(r— r') ( [l5| ) and Green functions ( [l7| ) provide the extension 
of the convenient multipole formalism towards sources with internal orbital structure. For freely propagating 
particle s \U {r) — 0], they naturally emerge in the partial wave decomposition of the scattering wave. (See 
Section nil.) 

Next, we turn our attention to the currents generated by multipole point sources. Assuming a superposi- 
tion of several of these sources at a fixed location r', a{r) — J2im ^imSim{^ — r')' '^he resulting scattering 
wave will read: 



(18) 



The current density j (r) due to this wave function then may be expressed as a bilinear form in the amplitudes 

Aim- 



Irn.l'm- 



(r)A, 



I'm' 



(19) 



Im I'm' 



where the vectors j 



lin.l'- 



' (r), given by: 



jirn.iw(r) {GL(r,r^£;)VGiw(r,r^£;)-Grm'(r,r^i?)VGL(r,r';£;)} , (20) 
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form the elements of the hermitian current density matrix. 

Of particular interest is the total current J{E) carried by the scattering wave ([l8|), i. e., the integrated 
current density j(r) (|l9|). Inserting the corresponding multipole source into the general current formula (Q), 
after integration by parts we first obtain: 



J{E) = - | lim 3 

ft r^r' 



E E KL(V)K,„,(V')G(r,r';i?) 



Im I'm' 



(21) 



Thus, total multipole currents can be extracted from the Green function G{r, r'; E) through differentiation 
with respect to r as well as r' and a subsequent limiting procedure, i. e., a number of elementary operations. 
We may further transform to obtain a more symmetrical representation of the current that again involves 
a hermitian matrix Jim,i'm'{E) [cf. (|l9|)]: 

— E E ^lmJlrn,l'm' {E)Xl'm' , (22) 
Im I'm' 

where the components of the total multipole current matrix Jim,i'm' {E) are given by: 

Jim,i'm'{E)^'- limKL(V)K,™(V'){G(r,r';i;)-G(r',r;£;)*} . (23) 

ft r^r' 

For simplicity, we will denote the (real and positive) diagonal elements of this matrix as the {I, m) multipole 
currents Jim{E): Jim{E) = Jimdm{E). 

If the Hamiltonian H possesses elements of angular symmetry, the number of non-vanishing total current 
matrix elements Jim.i'm'{E) ( p3| ) is greatly reduced. Assume that the generator of a rotation L commutes 
with H; then, it will also commute with the resolvent operator G — [E — H + \rj\^^ . Consequently, if the 
source states \(j) and \(j') are eigenstates of L with different eigenvalues, the mixed matrix element {a\G\(j'), 
and hence its contribution to the total current J{E) in (H), is bound to vanish. This implies that in case of 
full rotational symmetry of H, e. g. for a freely propagating wave with U (r) = 0, the total multipole current 
matrix (E3b is diagonal, Jim,i'm' (E) — Jim.{E)Sii'6„im', and the total current in (E2) becomes a simple sum: 



J{E) = J2hn {■^im\ Jim{E) (Scction 2.4). If the potential C/(r) is invariant merely with respect to rotations 
around the z-axis (like in the ballistic problem discussed in detail in Section orthogonality with respect 
to different values of m prevails: Jim.Vm'{E) = for m ^ to'. 

Finally, we set out to extend the Tersoff-Hamann description (H ) of the current to cover the case of 
multipole point sources. According to (Q), the current Jim{E) for a source located at R is formally given 
by: Jim{E) = 27r(R| K^„(p)(5(i? — Kim(p)|R)/?i. Expanding this expression again into a complete 
orthonormal set of eigenfunctions j^/^fi) of H, we find: 

Jlm{E) = ^E'^(^-^fl)|Kim(V)V'fi(R)|' • (24) 
fi 

Therefore, the multipole currents J/„j [E) ( ^ ) are proportional to the local density of the respective spherical 
tensor gradients of the eigenstates of the Hamiltonian H at energy E and position R. 



2.3 Analytical properties of multipole Green functions 

Our definition of multipole point sources (5/m(r — r') (|l5]) and corresponding Green functions Gim{j^, r'; E) 
( 1 7 ) by means of a differentiation procedure is a purely formal one. Thus, it appears appropriate to corroborate 
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their interpretation as sources of particles with definite angular momentum. Since we did not assume that the 
external potential V{r), and therefore the Hamiltonian H, commute with the angular momentum operators, 
the rotational symmetry of the multipole Green functions will generally be broken. However, in the immediate 
neighbourhood of the source, r — > r', the typical angular variation should manifest itself. 

Hence, in this section we study the local properties of Gjm(r, r'; E) in the vicinity of the point source and 
identify a universal behaviour that holds under quite general circumstances, viz., under the weak assumption 
that the potential V{r) is a locally analytic function. Unlike their global structure, the local characteristics 
of wave functions apparently have found little attention in the physics literature. The regularity properties of 
solutions to partial differential equations have been examined thoroughly, however, and the interested reader 
is referred to the monograph by Hormander [Hor63|. The subsequent results follow from an adaptation of 
this formalism to the multipole Schrodinger equation ( (1^ ) recently performed by one of the authors [ $ra99| ]. 

Mathematically spoken, the Hamiltonian operator H = T + V belongs to the class of elliptic operators 
[ Hor63 1, which implies that all solutions of the inhomogeneous problem [E ~ H]'ijj{r) = a{r) (||) are analytic 
in those sectors of space where the potential V{r) and the source term <T{r) are likewise analytic. It imme- 
diately follows that the difference of any two solutions to the point source problems in (|]) and (|l6|) will be a 
locally analytic function at r = r' (provided V{r) is regular at r'). Yet this statement does not apply to the 
Green functions themselves, since the source term is singular. To work around this problem, we employ the 

by a factorization approach 

We set G(r,r';i;) 



results of potential theory (Section 2.2) and eliminate the delta singularity in 
2MG(r, r')/(r, r'; E) /h^, where G(r, r') : 



l/47r|r — r'l is the Green function of 



potential theory. Then, /(r, r'; E) is a solution to: 



[R^Ar - 2R • Vr + R^k{Rf] /(r, r'; = 



(25) 



where R = r - r' and fc(R)2 = 2M[E - V{r)]/h'^. Additionally, f{r',r';E) = 1 must hold. In (||), 
the source location R = o is still a special (w eakly s ingular) point, but it can be proven that a unique 
solution exists that is locally analytic at r = r' | Bra99 1 . Furthermore, this function is real and symmetric 
[/(r, r'; E) — f{r', r; E)], and may be represented in the form (fc = fc(o)): 



/(r, r'; E) = cos (fc|r - r'|) + (r - r'yO{r, r') 



(26) 



(For free particles, the correction is absent. See Section 2.4.) Now, G(r, r'; E) is not the Green function itself, 
but differs from it only by an analytic expression ^(r, r'; E). Hence, G(r, r'; E) always may be displayed as 
a sum involving two regular functions: 



G{r,r';E) 



M f{r,r';E) 
2Tiff |r-r'| 



+ g{yy;E) 



(27) 



provided V{y) is analytic at r = r'. (Notably, (^) does not hold for the Green function of the Coulomb 
problem [Hos63, HP64|, where a logarithmic singularity is encountered.) Clearly, as r ^ r', the Green 
function (27|) adopts a universal shape regardless of V{r) and mimics G(r,r'): Locally, isotropy prevails. 
Obviously, the imaginary part in ( p7| ) remains well-behaved as r ^ r', thus justifying the expression ( pi] ) for 
the current J{E) generated by the point source. 

It is now but a small step to generalize these results towards the multipole point sources (5;„i(r — r') 



( 15 ). Applying the spherical tensor gradient to the Green function (E7h according to the rule ( 17 ), we obtain 



[Bra991: 



G, (rr'-E)- f^rn{r,r';E) 



(28) 
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where gim{^, r'; E) = Kim(V').g(r, r'; E) and fim{j^, r'; E) are again analytic functions at r = r', and the 
series expansion of the latter reads: 

hn{r,r';E) = Ki„,(R) [{kRy+^ni{kR)] +0{R'+^) . (29) 

Here, ni{z) denotes the irregular spherical Bessel function of order I [ A.S65[ |. In the limit r ^ r', the 
bracketed expression in (^) tends towards the constant value {21 — 1)!!, leading to the principal asymptotic 
form of Gim{j^, r'; E) in the vicinity of the source: 

GUry;E)^-^i2l^l)U^^^. (30) 

This expression, a multiple of the multipole potentials $/m(r, r') (|l2|) of potential theory, is again universal 
in scope. (Indeed, the relations ( pj| ) and (^4|) are special cases of the general result (^8|).) Obviously, the 
multipole Green functions (|l7[) are endowed with the desired property of local orbital symmetry. 



2.4 Example: Free multipole sources 



As an instructive example, we show how the mathematical developments of the previous sections operate for 
the simple case of freely propagating waves (C/(r) = 0). Here, the Green function of the scattering problem 
is just an outgoing spherical wave [ |Eco83| : 



M exp(ifc|r — r' 
27rft2 |r-r'| 



(31) 



with E = h^P/2M. We note here that G(f"'=)(r, r'; E) neady fits into the general scheme and (||), 
with analytical functions /(r, r'; E) ~ coskR, g{r, r'; E) = -M iiin{kR)/2iTh?R (w ith R — r — r'). For 
our purposes, a well-known expansion of ( ^T] ) into spherical Bessel functions is useful [ VIes64 1: 



ft2 



(32) 



where r< = min(r, r'), r> = max(r, r'), and j\{u), h^^\u) 
and Hankel functions of order A, respectively. 



n\ {u) + ij\ (li) denote the spherical Bessel 



With the help of (|32|), calculation of the multipole Green functions Gf^^'^\r , r' ; E) ( |17| ) and currents 
/m°°^ (i?) (^) becomes an easy task. We exploit translational symmetry and set r' = o. Then, we may 



replace jx{kr') ^ {kr'Y /{2l + 1)!! [ AS65|, and only the term {I, m) in the sum (M) survives: 



GfZ\B.^o,E) 



2Mk. 



(+) 



(fci?)r/„(i?) 



^^'^\\+\kR)YUR) 



(2; + !)!! 



K,„,(V')KL(r') 



27rft2 



(33) 



(The term in brackets is independent of m and yields fc'/47r.) Again, this expression conforms to the general 
structure of multipole Green functions laid out in ( |2^ and (^9|). From ( ^3|), we infer that G^^^°^ (r, r'; E) has 



the spatial symmetry of a spherical wave with (Z, m) orbital structure [ |Mes64 ]. Hence, the multipole source 
approach reproduces here the partial wave formalism of conventional scattering theory. 
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In a similar vein, we may extract the total multipole current jII^^^^ (E) ( p3| ) carried by a freely propagating 
wave. The action of the spherical tensor gradient on (BSh yields in the limit r ^ o: 



j; 



(free) 



iE) = 



4Mk 



{21 + 1) 



— K,„,(V')Kr,„(r') 



47r2ft3 



(34) 



The characteristic power law dependence of the near-threshold scattering rate on / is known as Wigner's law 
[ lWig48[ ]. 

It remains to identify the connection between realistic, extended sources cr(r) and the idealized point 
sources di„i{r — r') ( p3| ) introduced in Section 2.2, particularly for near-threshold scattering processes. 
Clearly, the replacement of aix) by multipole sources will only be feasible outside the source region, i. e., in 
the far-field limit, similar to the use of multipole potentials $;m(r, r') ^% in potential theory. Incidentally, 
any free-particle source a(v) can be substituted by a set of multipole sources (5;m(r — r') for an arbitrary 
choice of the source location r', as we are going to show now. With R = r — r', g(v) may be expanded into 
a spherical series at r': 



(7(r) =^Ki„(R)ai„(i?) 



(35) 



with even coefficient functions (Jim{R) accessible via projection over spheres centered in r'. For further 
reference, we also define related coefficents 7;m(r') by: 



,(r') = j dV Kr„, (r - r')a(r) = ARR^'+^aUR) ■ 



(36) 



Inserting the expressions (^2|) and (35 ) into the convolution integral (Q), and exploiting the orthonormaUty of 
the spherical harmonics, we obtain the multipole series for the scattering wave ^i^c'^'^\r): 



2Mk 



E 



di?' R"+^n{kR')ai^{R') 



,(+) 



ikR)YUR) , 



(37) 



(provided the source is entirely contained in the ball of radius R' < R). By comparison with the free 
multipole Green function (H), we find that V'sc = E;™ A,™(r'; E)Gf^^''\r, r'; E), where: 



\i^{v'-E) ^jrj, di?'i?"+2j,(fci?')'^;m(^') 



(38) 



Therefore, in the far field limit, the extended source a{r) may be replaced by a superposition of multipole 
sources Simi'r — r') ( p^ w ith strengths A/m(r'; E) (^8|). In the threshold regime (E ^ 0), we may replace 
the Bessel function in (38) by its limiting value, which causes the energy dependence to drop from these 
coefficients: 



\im{r';E = 0) = 



An 



(2/ + 1)!! 



7;™(r') 



(39) 



The source components then simply reflect the orbital structure of a{r) (^6|). (Although the components 
7im(r') generally depend on the choice of r', the leading non-vanishing coefficients of lowest order I are 
invariant with respect to shifts of the source location.) 
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Finally, we take a look at the total current J^^'"^'^\E) generated by the source cr(r). Starting from the 
general expression (Q), the expansions ( ^2| ) and ( ^5| ) yield the simple result: 

J^^'^'^iE) ^^Y.\^,„,ir';E)\V+' ^ \>^v^iy' ^ E)\\jfr\E) , (40) 

Ira Irn 

[cf. (0)]. In accordance with the discussion in Section 2^, the total multipole current matrix Ji^i'^m'i-E) 
(22), ( p3| ) for the free particle problem is diagonal, and J'-^'^°°\E) ( p^ ) is simply the properly weighted sum 
of the multipole currents (^4[). Concluding, we remark that near threshold, higher multipole contributions to 
j(ficc) (^^^ ^-g strongly suppressed. 



3 Ballistic multipole waves and currents 



In the preceding section, we applied the multipole wave approach to freely propagating particles. Obviously, 
the introduction of this new formalis m would be of little merit if it only served to reproduce the well-known 
results of partial wave decomposition [ Mes64 |. Therefore, we now tackle a less trivial problem, the propaga- 
tion of a multipole scattering wave accelerated by a homogeneous force field F = Fe^, i. e., in the presence 
of a linear potential U (r) — —Fz. (Subsequently, we will refer to the dynamics in this environment as quan- 
tum ballistic motion.) As an analytical expression for the corresponding ballistic Green function G'(r, r'; E) 
^ is available, the developments of Section 2J2 allow us to present the solution to the ballistic multipole 
problem in closed form. This section is devoted to a mathematical discussion of ballistic scattering, while the 
application of the multipole technique to physical phenomena under current study is deferred to Section ^|-^ 



3.1 Initial remarks 

As an extension of the method of Green functions G'(r, r'; E), the multipole source approach requires knowl- 
edge of the latter for a given background potential U{y). However, only for a scarce number of realistic, 
three-dimensional potentials the Green function is known in closed form. This selection evi dently include s 
free propagation (Section lA), but also c omprises the Coulomb potential ?7( r) = a/ | r — R| [|^os63 , |^P64 ], 
the isotropic harmonic oscil lator [^V71 ], th e homogeneous magnetic field [ |Gou72| , ^MM75| ], and parallel 
electric and magnetic fields [ Fab9 1 , KBKOl ]. For the b allistic problem, an anal ytic expression for the Green 
function was derived independently by several authors [ |DS76| , LLF90 , 3KK9 1 ] : 



G(r,r';£;) = 



M 



1 



[Ci(a+) Ai'(a_) - Ci'(a+) Ai(a_)] 



(41) 



where the arguments a± of the Airy functions Ai(u), C\{u) = Bi(it) + i Ai(it) [ |AS65t ] are given by: 

a± = -I3[2E + F{z + z') ± F\y - r' |] . (42) 



(A rather elementary derivation of this result is presented in [ BBG^98 1.) The parameter /3 in ( p^ denotes 
a characteristic inverse energy of the system (and /3F, correspondingly, an inverse length scale) that varies 
with the force strength F: 



P = {M/Ah^F^' 



,1/3 



(43) 
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Employing (3, it proves convenient for our calculations to introduce dimensionless quantities for the en- 
ergy, time, position, and momentum variables: 

p = PFv e = -2I3E 

K = k/(3F T = t/2h(3 (44) 

(The components of p will be denoted by p = (^j^jC)^- E. g., the arguments a± in ( |4l| ) read a± ~ 
e — (C + C') T \p ^ p'\ ) We use this notation in order to state two helpful integral representations of the 
Green function (ffll). The first follows from a Laplace transform of the time-dependent ballistic propagator 



[KBK02|: 



Gip, p';e)^ -2WFf / ,i(p-p')Vr+ir(c+C'-.)-i.Vi2 ^ (45) 

while the second form features the propagator in momentum space: 

d^K' e-'^'^-P+'^'-p') + k' + 2Te^) e''^'^-^'^^ '^^ . (46) 



Before starting out with the multipole calculations, we remark a general feature of the ballistic problem. 
From ( p4| ) and (p^, we infer that G(r, r'; E) is a functional of the variables r — r' and E + Fz' only. This 
translational symmetry is a consequence of the uniformity of the force field, where a shift of the source 
position merely alters the potential energy o f the emitted particles. We may take advantage of this invariance 
to relocate the source to the origin [ |KBK02| ]: 



G(y,y']E) = G{v-v'-o;E + Fz') . (47) 

Furthermore, the exclusive dependence on these variables allows us to replace derivatives of G'(r, r'; E) with 
respect to the source location r' by derivatives with respect to r and the energy E. One easily verifies that: 

d d d \ f d d d d\ 

We will summarizingly refer to the exchange of derivatives in (^8|) as the substitution rule for the ballistic 
Green function. 

3.2 Multipole Green functions 

We now characterize the multipole Green functions G;m(r, r'; E) ( [l7| ) of the ballistic problem. In view of 
the translational symmetry (p7|), we remark in advance that it suffices to discuss the case r' o; the general 
expression for Gim{r, r'; E) then follows by properly adjusting the positions r, r', and the particle energy E. 

By definition (p^), the multipole Green function is the spherical tensor gradient K/„i(V') of the ballistic 
Green function G(r, r' ; E) (p4|). Since this differential operator reduces to a polynomial in momentum space. 



it is advantageous to employ the integral representation (46) in performing the differentiation. In the hmit 
r' o, this procedure yields after a suitable shift in the remaining momentum integration variable q = 
K + TBz + 2p/t: 



Gi™(p, o; E) - dre'^'/^-'t'-^^^-^'/i^ J d'q Ki^ (q + re, - 2p/t) e' 



(49) 
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The latter integral involves the product of a Gaussian with a polynomial expression and hence allows for eval- 
uation in closed form. To this end, we first expand the shifted harmonic polynomial Kim (q + Tiz — 2p/r) 
into a spherical series with respect to q, as demonstrated in Appendix Q Since the Gauss ian p art is isotropic, 
the only contribution to the integral stems from the term with A = ^ = in the series ( |126| ), rendering the 
calculation trivial. Furthermore, by a repeated application of the translation theorem for harmonic polynomi- 
als (128), we may separate the spatial and temporal dependence in the argument of the remaining function 
Kimiriz ~ 2p/t). We then finally obtain a spherical series in p for the momentum integral in (S9b: 



d^q Kim (q + re^ - 2p/t) e" 



8^3/2 

(ir)3/: 



(50) 



where the translation coefficients Tjim are given by (129). At this point, it proves convenient to introduce a 
set of auxiliary functions Q^, (p, C; e) via: 



dr 



27r0F Jo (iT)fc+i/2 



exp < 1 



p2 



(51) 



Despite their rather involved appearance, for integer k these integrals can be systematically evaluated in 
closed form, yielding sums over products of Airy functions, as detailed in Appendix |B.1| . (Incidentally, 
apart from a constant p refactor, the function Q^jp, Ci e) equals the ballistic Green function in 2fc + 1 spatial 
dimensions [ |BBG+98| ]. E.g., a comparison with (p3 ) reveals that G{p,o;e) — — 4/3(/3_F)'^ Qj^(/9, C; e).) From 
(|49l)-(^l]), we then infer that the ballistic multipole Green function Gimiy, o; E) in the original coordinates 
is given by: 



Gim{v, o; E) = -4/3(/3F)'+3 ^ {2pFyTjim K,„(r) Q2^._;+i(/3Fr, /3Fz; -2(3E) . (52) 

j=\m\ 

The explicit expressions obtained from ( ^2[ ) quickly become cumbersome with increasing multipole order 
While the s-wave function G'oo(r, o; E) = G{r, o; E)/\/Att is displayed already in (p4|), here we merely 
state formulae for the p-waves in ballistic scattering [ Bra99[ ] : 



Gio{r,o;E) 



^ '^T- Ai'(a_) - Ci'(a+) Ai(a_)] 



+2pFr [P [z{2E + Fz) - Fr^] Ci(a+) Ai(a_) + z Ci'(a+) Ai'(a_)] } , (53) 



Gi,±i(r, o; E) = ^p^F^^^ { [Ci'(a+) Ai(a-) ~ Ci{a+) Ai'(a_)] 

-2pFr[Giia+)Ai{a-)+f3[2E + Fz]Ci{a+)Ai{a-)]} . (54) 

(The arguments a± of the Airy functions have been defined in (^2]).) The expressions for higher multipole 
order are, however, quickly calculated by means of (p2|). 

Despite being exact for all values of the parameters, the complicated structure of the explicit solutions 
( 52 ) to the ballistic multipole problem somewhat limits their practical use. Fortunately, comparatively simple 
asymptotic approximations for the vicinity of the source as well as for the far-field limit are available from 
(52). We begin with the source region, i. e., consider the case p ^ o. As explained in Section in this limit 
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the full multipole orbital symmetry of the corresponding source (5zm(r') (|l5]) should emerge in the ballistic 
multipole wave. This assertion is indeed true: Replacing the functions Q2j-i+i(p, C; e) in by their 
principal asymptotic form for p o (135), we conclude that the sum term with j ~ I becomes dominant in 
the neighbourhood of the source, and with the help of ( |43| ) and (129), we finally obtain the leading term in 

Gira{r,o;E) as r ^ o: 



K/„i(r) 



„2l + l 



(55) 



As expected, the ballistic multipole waves conform to the universally valid expansion (30) derived in Sec 



tion 2.3 



In the far-field region z oo, the ballistic multipole waves assume a fairly different character. Since 
the presence of the uniform force field manifestly breaks the full rotational symmetry incorporated in ([T^, 
the characteristic angular dependence of the scattering waves propagating in a central potential U{r), as 
exposed e. g. by the free spherical partial waves (p3|), is destroyed. (Note, however, that the uniform field 
problem still retains an axis of rotational symmetry. Consequently, the ballistic Green functions G;,,, (r, o; E) 
(52) are eigenfunctions of the angular momentum operator L^, and terms with different values of m don't 
mix.) For large distances from the source, the external potential bends the scattering waves in the direction 
of force, and the particle current forms a beam centered around the z-axis that continues to spread in the 
lateral directions while its density profile settles into an invariant shape that reflects the orbital structure of 
the source. In a sense, therefore, the force field projects the spatial particle distribution near the source 
onto an enlarged distant plane image. (For electrons in a homogeneous force field, this effect has been 
dubbed "photodetachment microscopy" by Blondel et al. [BDD96, BDDV99|.) We are thus interested in the 
asymptotic shape of the ballistic multipole wave in the vicinity of the z-axis, i. e., in the limit ( ^ oo while 
p — ( remains finite. (The exact formula (52) is of limited use in this region due to large-scale cancellation 



between the various contributions to the sum.) Fortunately, there exists a systemat ic appr oach (employing a 
saddle-point approximation) to evaluate all ballistic waves in the far-field sector | Bra99 1 . We refrain from 
displaying it here; instead, we only cite its rather simple result that involves again a harmonic polynomial 
operator: 



Glrnir,0]E) 



2 



K 



Im 



(3Fx (3Fy . d 



9a _ 



Ai(a_) 



(56) 



(Here, the argument a+ ( p^ approaches —oo, while q:_ remains bound.) Note that all multipole 

waves share a common propagating factor, whereas the asymptotic beam profile is contained in the derivative 
of Ai(a_). Explicit expressions for s- and p-waves are tabulated in Section Q 



3.3 Ballistic multipole currents 

Despite being a fundamental set of quantities, the ballistic multipole Green functions G/m(r, o;i?) 
are not directly accessible in experiment. This privilege is actually reserved for the multipole currents 
j/m,;'m' (r, o; E) (|^ and Jhn,i'm' (E) (f3|) which we discuss next. 

We start out with a brief analysis of the current density j(r) ( |l9|) a ssociated to a ballistic scattering wave. 
This quantity is accessible from the Green functions G/m(r, o; E) (|52[) through differentiation, and the calcu- 
lation of the exact matrix elements j/m,;'m' (r, o;E) is tedious yet straightforward. The resulting expressions 
are lengthy, and we will refrain from displaying them here except to notice that the current density compo- 



nent Joo'oo(r, o; E) in field direction for a ballistic s-wave has been published before [tCBK02|. However, 
the current density profile in the far-field section has been assessed experimentally [ BDD96| , ^DDV99 ], so 
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it appears worthwhile to deduce asymptotic expressions for the various multipole current matrix elements 
in this limit. Since the scattering wave in this regime closely follows the direction of force, we will only 
examine the z component of the current density vector, i. e., confine our attention to (r, o; £'). It 

has been argued [ |BDDV9^ ] that for z oo, the velocity of the particles in the matter wave is determined 
by their uniform acceleration in the field, regardless of the details of the emission process, and hence is 
given by v{z) ^ ^ 2Fz /M. In this approximate picture, the resulting current distribution j'-^^(r, o; E) is 
simply proportional to the particle density p{v) = |-0sc(r)P of the scattering wave: ^''•^■'(r) ^ p{v)v{z). 
This notion is indeed corroborated by a more accurate inspection [ |3ra99 |: Using the Wronskian relation 
Ci(z)* Ci'(z) - Ci{z) Ci'(z)* = -2i/7r [ |AS65| , one quickly arrives from (^6|) at the far-field current distri- 
bution (z oo): 



E) 



K 



Ira 



(3Fx (3Fy 



(2/3F)'+''+^i'' 



^— ' dc 



Ai(a_) 



f3Fx PFy . d 



^—a^ ^— ' dc 



Ai(a_] 



(57) 



Here, «_ ( |4^ is bound, while a+ —oo. The resulting expressions for ^ = 0, 1 are utilized in Section 4.2. 

In the second step, we determine the total current J{E) ( pi] ) carried by a ballistic scattering wave. Its 
matrix elements Jim.i'm' [E) ( ^3| ) follow from the G reen function G(r, r'; E) (^Tj) by differentiation and a 
subsequent limiting process, as discussed in Section 2.2. Before heading towards the calculation of those 
currents, we note that due to the remaining rotational symmetry of the force field, the total current matrix is 
diagonal with respect to the quantum number to, i. e., Jimj'm' (E) = for m to'. (See Section 2.3.) 

As the multipole formalism invokes spherical tensor gradients K/„j(V) (^3|), it is again favorable to 
perform the actual calculation in the momentum-space representation of the Green function G(r, r'; E) (|4^. 
The differentiation and limit operations then become trivial, and the matrix elements of the ballistic total 
current read in integral form: 



l-l' + l 



(58) 



In the next step, the momentum and temporal contributions in the arguments of the harmonic polynomials 
appearing here are disentangled by means of the translation theorem (128) (see Appendix This renders 
the angular integration in ( p8| ) straightforward, and the remaining momentum integral is of Gaussian type, 
leaving only the temporal integral to be evaluated. It may be displayed as a sum over the auxiliary functions 
Qikie) closely related to the expression appearing in the preceding section: 



Qifc(e)= lim 3[Q,(p,C;e)] 



dr 



-1/2 



-iT712 



(59) 



(A similar set of functions is discussed in [ [VlFSFOd[ |.) Their properties, in particular, their resolution into 
products of Airy functions, are discussed in Appendix B.2. With the definition (]59|), the total ballistic multi- 
pole currents Jim,i'rn'{E) (ESb finally read: 



Hm,l'm' 



(E) 



M 



min(/,^ ) 
j = \m\ 
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For ^ = 0, 1, these currents are listed in Section 4.3 



Despite being an exact expression, it is not easy to comprehend the behavior of the total multipole currents 
from (pOl). To elucidate it, we proceed with an asymptotic expansion of the multipole currents Ji„i{E) (here, 
we treat the diagonal matrix elements only) in terms of the energy parameter e ~ —2(3E valid in the sector 
|e| 3> 1. Clearly, the problem reduces to finding asymptotic forms for the auxiliary functions Qi/i.(e) (|5^, 
and its somewhat intriguing solution is displayed in Appendix B_3. Here, we merely state the subsequent 
results for the ballistic current Ji,n{E). 

We have to distinguish between two mathematically and physically rather different cases. To start out, 
we consider the limit e ^ 1 pertaining to large negative energies. This corresponds to the interesting phe- 
nomenon of ballistic tunneling from a point source that has no classical counterpart. The leading asymptotic 
term for the total multipole current Ji„i{E) then reads: 



Jlm{,E) 



Mk2'+i (2; + l)(; + |m|)! 
47r2ft3 \m\\{l ~ \m\)\ 



I3F 

K 



3|m|+3 



-K^/6(/3F)^ 



(61) 



Here, k = 2fiF^ denotes the evanescent particle momentum at the source. As expected, tunneling is 
exponentially suppressed. We also note that for fixed I, the current strength declines with increasing quantum 
number |r7i|. This is evidence for the centrifugal suppression of tunneling discussed below. 

In the opposite limit of classical ballistic motion, k oo (where k = 2f3Fy/^, e ^ —1), the multipole 
current Jim{E) may be interpreted as a sum of two terms of different character: The dominant secular 
contribution is independent of the field strength F and thus obviously must reproduce the total current for 
freely propagating partial waves (Section 2.4) given by Wigner's law (p4|). The free-particle expression is 
then modified by an oscillating contribution akin to (|6l|): 



Jlm{E) 



Mk 



21 + 1 



47r2ft3 



,2(2? + l)(/ + |m|)! 
\m\\{l ~ \m\)\ 



k 



3|m|+3 



COS 



1 

6 \I3F 



(62) 



The modulation becomes most apparent for linear polarization (m = 0), but even then it is suppressed with 
respect to the secular part by a factor k^^. However, due to its rapidly oscillating behavior it still may imprint 
a conspicuous pattern onto the general trend of Ji„i (E): For m = 0, the energy derivative of the current reads 
asymptotically: 



2L2i-l 



(2^+1) 



27r2n5 



■ sm 



|(2/3i?)3/^ + (-l)'J 



1 2 



(63) 



Thus, the current J;o {E) grows monotonically with E, yet becomes stationary at a series of energy values 
Ei,i — [3tt{4i> + 21 — l)]2/3/8/3 {v integer, v > —1/2), rendering a "staircase" structure in plots of Jiq{E) 
versus the energy E. 



3.4 Semiclassical theory 

In the previous sections, we derived the exact quantum mechanical solution to the ballistic multipole prob- 
lem as stated in equations (|2|), (^7|), and (^0[), respectively. Although a major feature of our presentation, 
the physical content of these mathematical expressions is not extracted easily. Hence, in the following para- 
graphs, we will provide some surprisingly simple arguments of essentially semiclassical nature that within 
their range of vaUdity faithfully reproduce the features of the quantum solution. As before, we start out with 
the current distribution in the far-field sector (z oo). 
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Like in the preceding section, we have to distinguish two radically different energy regimes. In quantum 
motion, the initial particle energy E can be negative, leading to ballistic tunneling, which is a phenomenon 
with no immediate classical counterpart. We will assess it below. For the moment, however, we consider 
the simpler case of classically allowed ballistic motion originatin g from a point source, a problem alread y 
studied by Galilei. Our approach, like some earlier contiibutions [ [FabSlj |DK082| , |Du89i |KO90b| , |BDD96| ], 
is bui lt upon th e Hamilton-Jacobi theory of uniformly accelerated motion which is discussed in detail in 
Ref. [ BBG+98 1. We assume that classical bodies of some fixed energy E are emitted at the origin r' = o and 
travel towards a destination position r on a distant screen perpendicular to the direction of force (z — const., 
z oo). The freely falling bodies then may follow two different parabolic tracks: The "fast" trajectory, 
denoted by (— ), directly leads to the location r, whereas the "reflected" or "slow" trajectory (+) first grazes 
the turning surface, a rotational paraboloid with the source as its focal point given by a_ = (||) which 



presents the maximum range of classically allowed motion [BBG+98], before arriving at r. (The situation is 
depicted in Figure [l|.) On the screen the particles will populate the circular disc a- < whose radius Rc\ is 

AEz/ F. Correspondingly, as z ^ oo. 



asymptotically given by 



a- - -2(3E (1 - R^/Rl 



(64) 



Thus, for large distances the current profile settles into an invariant shape of absolute size i?ci- 

Next, we discuss the projection properties of the force field for large z. The relative lateral position 
R/Rc\ on the screen is determined by the angle of emission 9 relative to the direction of force. Since motion 
perpendicular to F is uniform, R is the product of the corresponding component of the initial particle velocity 
with the time of flight: R — T{9; E)viSm6. However, in the far-field limit T{9;E) is asymptotically 
independent of 9 and E, and we find the simple relationship R{9) = Rd sin 9. Hence, the cone 6 ~ const, will 
be projected onto a circle of radius Rd sin 9, and trajectories starting under opposite angles {9, (f), (tt — 6*, 0) 
will asymptotically share the same destination r on the screen. 

The distribution of emission angles at the source and the current profile recorded at the detecting screen 
are linked through the classical differential cross section dad/dVl of the ballistic problem in the asymptotic 
limit. For both the direct and reflected path respectively, we find: 



dod 

on 



{R;E) 



RdR 
im9d9 



Rd cos 9 



R'dJl-RyU; 



(65) 



The den sity of trajectories dfl/d(Td{R;E) diverges at Rd, as expected near the caustic surface a- ~ 
[ Sch81 1. In the actual current density distribution, this divergence will be smoothened out into a bright 
fringe. (If, however, emission perpendicular to F is suppressed, as e. g. for a multipole source, the edge 
of the current spot appears quenched.) 

In an entirely classical picture, the differential cross section ( |6^ ) would completely determine the particle 
distribution at large z. However, a fundamental modfication is imposed by the wavelike nature of quanta. A 
freely falling particle of fixed energy E will travel along both the "slow" and "fast" paths, but will accumulate 
different phases in the process, leading to interference between the particle waves. As the phase difference 
depends on the location r — at R — Rd the classical trajectories will coincide, whereas maximum phase 
difference occurs for particles emitted (anti-)parallel to the direction of force — constructive and destructive 
interference will alternately take place, leading to a circular interference pattern superimposed onto the cross 
section (pS]). Conceptionally, the situation resembles the traditional double-slit setup |MJ59, FH65| used to 
discuss two-path interference. From a practical point of view, the ballistic problem enjoys the advantage 
that no mechanical "slit" is involved as the selection of paths is solely accomplishe d by the force field . The 
predicted interference patterns were first observed experimentally by Blondel et al. [ BDD96 , ^DDV99| . 

A quantitative semiclassical analysis rests on the quantum phases cr± {R, (j); E) carried by the interfering 
trajectories (+), (— ). Generally, three different terms contribute to (7±{R, cj); E) (see also Figure |l]): First, 
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Figure 1: Semiclassical ballistic motion in the far-field limit. Within the disc R < Rd, two parabolic paths 
(bold) emitted under equal angles d,TT — 9 from the multipole source will join the origin with the destination. 
Unlike the fast path (+), the slow path (-) undergoes a reflection at the parabolic turning surface a- — 0. The 
accumulated phases, together with their initial "atomic" phases inherited from the point source, determine 
the exact shape of the interference pattern on the screen. 



the trajectories will "inherit" atomic phases ^{9, (p), 7(7r — 9, cj)) deriving from the angular phase distribu- 
tion of the point source; second, particles traveling along the classical paths will gather dynamical phases 
which semiclassically are determined by the reduced action W^^\r, o; E) along the respective trajectories 
||BBG+98| which reads in terms of a± (p^): 



(±) 



2^ 



)3/2±(-a.) 



3/2 



(66) 



Finally, the "slow" trajectory (+) undergoes an additional "phase jump" of — 7r/2 due to reflection at the 
turning surface [Sch81|. Hence, the semiclassical phases (t±{R, (p; E) along the trajectories (+), (— ) read: 

1 



<7+{R,4>-E) = 
a^{R,4>;E) 



1{ 



\w^^{\v,o-E). 



(67) 
(68) 



We now put together the semiclassical theory of ballistic motion in the limit z ^ oo. Denoting the angular 
amplitude distribution at the source by A{9, (p) — \A{9, (p) \ exp[i7(6', (p)], the semiclassical current density 
distribution j't'^^ (R, (p; E) on the screen will be given by the differential cross section d<7ci/dfl{R; E) 
accounting for the projection properties of ballistic motion, modulated by an oscillating term representing the 
combined effects of source emission characteristics and quantum interference (where sin 9 = R/Rc\): 



jt^)(i?, 0; E) = ^{R; E) \A{i, - 9, 0)|e-+(^'^^^) + </>)|e-^(^^*^^) 

0(Jc\ 



(69) 
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In general, little information about the emission characteristics can be extracted from the current profile. 
However, the contributions of source structure and quantum interference will disentangle if the emission 
pattern shows reflection symmetry with regard to the x — y plane, i. e., \A{tt — 9, 0)p = \A{6, (f))\'^ holds. 
I reduces to a product of three independent factors [ ^K082 |: 



Then, 



ji^'\R,cb;E)^4\A{e,cf>)[' 



dac] 



{R; E) cos^ 



\R,(t)-E) - (7_{R, 0; E) 



(70) 



This condition always applies for pure multipole sources (5;m(r) (Uj). We now calculate the semiclassi- 



cal approximation i'f^\^ 



(r, o; E) to the ensuing multipole current density distribution im(^^ °' ^) (|57|). 



Apart from an overall factor Ao{E) the angular amplitude distribution of this source is given by a spherical 
harmonic function: A{9,4') = Ao{E)Yi„i{d,(j)). In the classical picture, the total emission rate J;„j(£^) 
(60), and thus Ao{E), remains undefined. Here, we pragmatically adopt the value quantum theory pro- 
vides for field-free emission, i. e., the Wigner law (|34 ) from Section 2.4, and set Aa{E)'^ = jj^^'^^^E) = 
^y^^2;+i^^^2^3 Inserting an explicit formula [ [Edm57 | for Yim{6, (f)) into ([70|), eliminating the angle 6 and 
approximating the semiclassical phase a±{R, (t>',E) (|66|)-(p8|) using (pSb we obtain after a short calculation: 



i£L(r,o;£;) 




(71) 



Here, P"^{z) denotes the associated Legendre polynomial [AS65|. Multipole sources thus generate two 
different patterns of sharp ringlike interference fringes. The upper sign in the interference term applies for 
even parity, i. e., for even I — |m|, whereas the lower sign is valid for odd I — \ni\. Therefore, under change 
of parity, the interference ring pattern will reverse. 

With increasing energy E, new interference fringes spring up in the center of the current distribution. 
There, the phase difference (7+ — cr_ ( |67[ ) is linked to the action functional along the "(+)" path segment 
returning to the source, known as the "closed orbit": WcoiE) = 4/i(2/3i?)'^/^/3. It has been noted |DD891 
that destructive interference in the closed orbit takes place whenever E = Ei,i (see Section 3.3). Thus, the 
stationary points in the total current Jiq{E) ( |63| ) coincide with the appearance of a new dark interference 
fringe in the pattern, providing a crude explanation for the staircase structure prominent in the multipole 
currents with m — Q. 

We now turn to the case of ballistic tunneling (E < 0). It is of considerable theoretical interest in being 
probably the only example of a tunneling process in three spatial dimensions yielding to an analytic so lution. 



and thus forms an importa nt testing ground for semiclassical theories of multi-dimensional tunneling [ |KP37 
HS67 , BBW73 , HFCK90 1 . Unlike their one-dimensional textbook counterpart, realistic tunneling problems 
are notoriously intricate and show a wealth of surprising features (e. g., trajectories in co mplex spac e) absent 
in the former For a discussion of ,s-wave ballistic tunneling, we refer the reader to Ref. [BBG+98]. 



A systematic study of ballistic multipole tunneling is outside the scope of this article [ Bra99 1 . Rather, we 
present a simple heuristic approach that might be dubbed "analytic continuation": Even as the notion of a 
trajectory becomes dubious when applied to tunneling, we formally take over the results (p3)-(71 ) achieved 



above for classically allowed motion, subjecting them only to minor modifications. We note that the square of 
the classical radius as defined by a_ 0, is now formally negative. We replace it by its absolute value 
= 4\E\z/ F. Similarly, quantities derived from it will change accordingly; in particular, we 
~3), dcjt^n/d^{R, E) = + ^V-R?u„ (H) for the "cross section," 



p2 _ 
■"'tun — 

find q;_ 



-2/3i?(l + i?V<J: 
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and the "projection law" cos 9 = -^/l + R? / R^^^. Clearly, these expressions are devoid of any geometrical 
interpreta tion. More over, the tunneling action functional, and hence the phase along the "trajectory," becomes 
complex | ]BBG+98t ]: 



Wtu„(r,o;i;) 



2h 



(72) 



Since the tunneling wave function must decay exponentially, only one of the two possible extensions of ( 
is physically acceptable. This restriction lifts the path ambiguity, so interference is absent in the tunneling 
regime. Therefore, (|69|) is replaced by: 



CCTtun V O 



(73) 



For the special case of a pure multipole source 5im (r) (|15|), we finally obtain with the evanescent momentum 
Hk = yj2M\El in analogy to (0: 



.ItunI , „^ 

jL,;m(i-,o;£;) 



Mk2'+i 21 + 1 {I -my. 




exp 



-2(3E f ] 



3/2" 



(74) 



(For real z > 1, the Legendre polynomial P^iz) is generally complex and double-valued. The modulus 
|P™(z)p remains unaffected by this complication, however.) 

It is instructive to examine the paraxial limit i? of (Q). This yields the following approximation to 



JiTiL (r, o; E), vahd for small lateral distances R [ |BBG+98| 



Mk 



.21+1 



(2/ + !)(/+ |m|)! 



\m\)\ n2|m|+2 
' ^tun 



exp 



kR^ 

' 2z 



(75) 



This formula allows for a simple interpretation as a product. The exponential term exp[— K'^/6(/3i^)^] equals 
the WKB penetration factor for a one-dimensional linear potential ramp, while the prefactor (i?/-Rtun)^'™' 
covers the effect of centrifugal repulsion. It is worth noting that the tunneling current distribution, apart 
from this prefactor, possesses approximately Gaussian form: J'"" (i?, 0; £J) oc exp(— Ki?'^/2z): Remark- 
ably, unlike the total tunneling current Jim{E), the shape of the lateral current profile is independe nt of th e 
force strength F. This prediction is experimentally confirmed in field emission from ultrasharp tips [ HM93 ]. 
We conclude our investigation by pointing out that integration over the current profile (^5|) provides an in- 
dependent means to calculate the total ballistic multipole current Ji„i (E) in the tunneling limit. The result 
coincides with our earlier finding (|6l|). 



4 Application to photodetachment processes 

The multipole wave functions and currents characterizing quantum balUstic motion transcend merely math- 
ematical interest. Remarkably, these quantities have recently been assessed experimentally in a rather direct 
fashion, studying electrons accelerated in a homogeneous electric field and ultracold atoms subject to the 
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gravitational force, respectively. In this section, we start out with a short description of near-threshold pho- 
todetachment of negative ions in a the presence of an external uniform field. This problem has been studied 
extens i vely theor etically (see [ |Fab8ll |Fab82i |DK082i [WRGSSl |DD88i |DD89i ^3u89l |KO90"bi |K090H, |Fab94 , 
Gol97|, |MFSFOO|| and references therein) as well as experimentally [ |BMS+87| , |GDL93b| , |GDL93a| , |BDD96 , 
BDDV99, GGM+01]. Being aware that many of the subsequent results have appeared in the literature, we 



nevertheless chose to discuss the multipole source model of photodetachment as it presents the most coherent, 
and by far the simplest, description of the effect. However, we strive to emphasize the less known and probed 
aspects of the subject. For a brief introduction to the source formalism for s-wave photodetachment in an 
electric field environment, we refer the reader to our previous article [ KBK02 ] . 



4.1 Photodetachment as a source problem 

In the photodetachment setup, a beam of negatively charged ions traverses the focus of a laser beam whose 
frequency uj closely matches the electron affinity of the ion, i. e. the binding energy Eq of the excess charge. 
Some ions absorb a laser photon and subsequently emit an electron of energy E = huj — Eq into a con- 
tinuum state to become a neutral atom. In the presence of a homogeneous electric field F = — eE, these 
electrons are accelerated towards either a counter, allowing measurement of the total photocurrent J{E) 
[GDL93b, GDL93a GGM+Ol"| |, or a spatially resolving detector plate that records an image of the photo- 



electron distribution jz (r, o; E) [BDD96, BDDV99|. (For technical reasons, the ion and laser beams and the 



electric field are usually oriented mutually perpendicular We choose the z-axis as direction of force.) 

In the source formalism (Section 2.1), the photodetachment phenomenon is interpreted as the scattering 
of the ionic electrons at the quantized electromagnetic laser field A, i. e. the coupling term W = — e?i^(p • 
A)/Mc is treated as the interaction potential. Assuming that the external electric field and the laser field are 
sufficiently weak not to disturb the electronic configurations of the initial ionic ground state jV'ion) and the 
emerging neutral atom in its ground state iV'atom), we proceed to project these states, leaving an effective 
inhomogeneous Schrodinger equation for the detached electron (0): 



E 



(76) 



where C/atoin(i") represents the short-range interaction between the emitted electron and the remaining neu- 
tral atom. In leading approximation, this potential can be neglected. (In the related photoionization effect 



I KO90b , NOV+021, a long-range Coulomb attraction between the emitted electron and the emerging ion pre- 



vails which must be included in the external potential: Uo^t (r) = Fz + e /r. The additional term renders the 
treatment of near-threshold photoionization considerably more difficult.) In strong electric fields, the effects 
of J7atom(r) may become important and must be included in a perturbative rescattering series. Similarly, the 
influence of the oscillating laser field on the electronic motion may be taken into account by a Floquet series 
expansion. (The electron dynamics in the laser field becomes dominant only for extreme intensities where it 
causes a wealth of new phenomena, e. g. high-harmonic generation [Gav92 LKKB97].) A comprehensive 



study of these corrections is performed in Ref. [MFSFOOj. However, for moderate fields and laser intensities 



and frequencies close to the detachment threshold, they may be safely ignored. 

In order to reduce the problem to fit the ballistic multipole formalism discussed in Section ^ we finally 
approximate the source term a{r) — ('0atoin|W^|V'ion) in (76) by a properly chosen multipole point source 



(15). To this end, we first note that a(r) is limited in extension to the size of the parent ion and thus is much 
smaller than the initial wavelength of the photoelectron. Similarly to the case of free-particle sources (Sec- 
tion 2A ), the details of the source structure then may be condensed into a single parameter, the source strength 
C. (Otherwise, finite-size effects have to be taken into consideration. For Gaussian sources, the deviation 
from the point source behaviour has been studied in detail [KBK02]. See also Sections.) This leaves only 
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the multipole character of the source to be determined, which in turn is fixed by the selection rules for dipole 
radiation. Since the emission into channels of higher angular momentum at energies close to threshold is 
strongly suppressed (as exemplified e. g. by the Wigner law (Q) valid for freely propagating particles), only 
the lowest permissible multipole order is appreciably popula ted. In most cases, t his effect leads to i sotropic 
emission of the photoelectron (I = 0) from the point source [ |BDD96| , |BDDV99| , |GGM+01| , |KBK02| ]. How- 
ever, if both the parent ion and the emerging neutral atom possess S ground states, the photoelectron must 
carry the spin of the absorbed laser photon and is therefore emitted into a p-wave, where the distribution onto 
the various m sublevels is determined by the laser polarization vector e: 



E 



2M 



Fz 



^(r) = C(e. V)<5(r) ^ Ai,„,Si„(r) 



(77) 



(The most prom inent member of this class is the hydrogen ion H~ first studied experimentally by Bryant 
et al. [BMS_J_87|.) In this form, the photodetachment problem immediately yields to a description in terms of 
the ballistic multipole waves Gio(r,o;£') ( ^s] ) and G'i_±i(r, o; (^). Some results for the ensuing current 
distributions are presented below. 



4.2 The far-field current profile 



According to the semiclassical model of quantum ballistic motion (Section 3.4), the uniformly accelerated 



wave function is characterized by a circular fringe pattern imprinted on the Green function by two-path 
interference in the force field (Figu re p]), an d approximations similar to were first published by Fabrikant 



|Fab81, pab82 |. Demkov et al. [DK082] not only improved on these results, but also realized that the 



interference pattern should be experi mentally observable in near-threshold photodetachment. This assertion 
was finally verified by Blondel et al |BDD96, BDDV99 I in a seminal series of experiments: They recorded 
the photoelectron distribution generated by various ion beams in a homogeneous electric field with a spatially 
sensitive detector plate and in the course established a new method for the precise determination of electron 
affinities ("photodetachment microscopy"). The observed interference patterns are truly macroscopic; for 
electrical fields E of a few hundred V/m and a source-detector distance z = 0.514 m, the fringe diameters 
are in the mm range. However, all experiments were performed on ions that detach electrons into s-waves. 
Due to the low photoabsorption rates imposed by the Wigner law (^4|), photocurrent profiles in near-threshold 
p-wave detachment are more difficult to measure and so far have not been recorded. However, source theory 
predicts an interesting dependence ( [tt] ) of these patterns on the laser polarization e that is naturally lacking 
in isotropic emission. We discuss some examples below. 

First, however, we briefly comment on s-wave photodetachment. In this case, the electron is effectively 
emitted by an isotropic point source a{r) = C5oo(r), and its wave function is thus proportional to the 
ballistic Green function G(r, o; E) (^l]). The density of electrons collected on the distant detector hence 
serves as a measure of the current density j^^^ (r, o; E) assigned to the Green function (^ which is stated in 
exact form in Ref. |KBK02|, and is in very good agreement with experiment. Yet, due to the considerable 
distance of source and detector, a simplified analysis using the far-field approximations (|6|) and ( ^7| ) for the 
ballistic Green function and its current density is equally suited. For s-wave emission from a unit strength 
source (C = 1), we obtain: 



Goo{r,o;E)^Ai(3{l3Ff 



Ai(a_) , 



(78) 
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(where a± is defined in (|4^), while the ensuing current density matrix element reads: 

j^l],„{r, o; E) ^ Ai(a_)2 , (79) 

Thus, for fixed distance z, the electron density varies with the radial distance R like the square of the Airy 
function Ai[-2p E{l - / R%)f (|3), where Rl^ = AEz/F (see also Section K^. This formula was used 



by Blondel et al. [ BDDV99 1 in the analysis of their photodetachment experiment. 

Now we turn to the case of p-wave photodetachment. From ([77|), we expect that the electron wave func- 
tion i^{y) is a superposition of the I = 1 multipole Green functions Gio(r, o; E) ( ^3] ) and G'i,±i(r, o; E) ( p4[ ) 
whose relative weights depend on the laser polarization e. For our purposes, again their far-field asymptotics 



(56) suffice: 



Gio(r, o; E) ^ -8\/3i/3(/3F)4-^l^±L Ai'(a_) , (80) 
Gi,±i(r, o; E) ^ T4^/3/3(/3J-)^ (x ± iy) Ai(a_) . (81) 



Their associated current density matrix elements jj^j^ j^^^, (r, o; E) are easily evaluated from eq. (|5^). 

Here, we examine four specific setups of the photodetachment experiment in detail. As noted before, 
the ion and laser beams, and the direction of force should be in mutually orthogonal orientation; we choose 
them as the x-, y-, and z-directions of our coordinate system, respectively. Then, the polarization vector 
e is confined to the x — z-plane. Some natural choices for it are e^r = (parallel or 7r-polarization with 
respect to the electric field E), — Bx (perpendicular or a-orientation), Ccirc = (e-z + iea:)/-\/2 (circular 
polarization), and ctiit ~ (e^ + ex)/V2 (linearly polarized under an angle of 45° to E). The corresponding 
multipole sources in ([77|) then read, apart from a constant of proportionality: 



cr,r(r) = (5io(r) , CTcircW = i[-i(5ii(r) + \/2(5io(r) + i5i,-i(r)] 

1 

V2' 



(7,(r) = ^[-^iiW - <5i,-i(r)] , (Ttiit(r) = ^[-Sn{r) + V2(5io(r) + <5i,_i(r)] . (82) 



Note that crcirc(r) and crtiit(r) only differ in the relative phase of their multipole components. A brief calcu- 
lation using ( |57| ) then yields the asymptotic photocurrent density profiles generated by these p-wave sources: 

ji^Hr,o;E) - ltf^\ Ai'{a.r , (83) 
7r^n,(— a+j 



jW(r,o;i?) ^ Ai(a_)2 , (84) 



ji-l{r,o;E) ^ llf^' fAi'(a_) - ^ Mia^))' . (85) 



j[:^{r,o;E)^ fAi'(a_)^ + i^Ai(a_)^) , (86) 
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(In a differe nt, rath er complicated manner, approximations equivalent to (83) and (|^ were derived by 
Golovinskii [ |Go197| ].) 

We now highlight some properties of these d istributio ns. They are depicted in Figure ^ for the set of 
parameters actually used in Blondel's experiment [ BDD96 |, viz., initial electronic energy E = 6.08- 10~^ eV, 
electric field strength E = 116 V/m, and detector distance z = 0.514 m. A concentric arrangement of the 
interference rings only occurs for tt- and cr-polarization of the laser beam. However, the nodes and maxima 
are interchanged in the respective images. This effect was already predicted in Section 3.4 and may be 
considered as a hallmark of the differing atomic phase 7(0, 0) of the interfering trajectories. The fact that 
the fringe pattern encodes information about the orbital structure of the source lead Demkov et al. [ pK082| | 
to its interpretation as an image of the atomic wavefunction, generated by "photodetachment microscopy." 



(In (T-polarization, the underlying interference pattern is the same as in s-wave photodetachment [ BDDV99 



KBK02].) Likewise, a sharply defined interference structure is observed for circular polarization ech 
as this source conforms to the condition outlined in (^0|). However, the mirror symmetry x — > —x, present in 
the other three plots, is conspicuously broken, as the number of fringes differs on the left and right side of the 
current density profile. Finally, due to differing weights \A{9, (p)\ of the interfering trajectories, no clear-cut 
ring pattern is present in the tilted case ctiit. The corresponding image rather looks like a combination of the 
patterns obtained in tt- and cr-polarization. Indeed, eq. (p3l)-(|8^ show that j^^-^l (r, o; E) = [ji^'' (r, o; E) + 
ji^\r,o;E)]/2ho\ds. 



4.3 The total photocurrent 

For historical reasons, the total detachment rate J{E) of photoelectrons as a function of the electron excess 
energy E attracted more attention than the photocurrent distribution j(^)(r, o; E), despite the fact that J{E) is 
available from the latter quantity by integration. Stirred by the experiment of Bryant et al. | BMS+87 1 , a series 
of theoretical articles approached the problem, generally based either on integral representations of the Green 
function (0) [ |WRG88[ , |Fab94| ] or an analysis using Fermi's golden rule (|T^ | pD88| , pD89( ]. The theory of 



ballistic multipole source developed in this article yields a closed expression for all multipole photocurrents 
Jim{E) (|6C|), thus obliterating the need for the rather involved calculations in the former approaches. 

We again briefly comment on the cas e of isot ropic emission of the photocurrent. For a source of unit 
strength (C = 1), we immediately obtain [ KBK02 1 : 



27rri3 



eAi(6) 



(87) 



where e = ~2(3E (|3). This result was also fo und by F abrikant [ [Fab94| , and, via the density of states rep- 
resentation (|ll]), by Luc-Koenig and Bachelier [ LKB80 |. The near-threshold s-wave photocurren t spectrum 
has been recorded with high precision in the experiments by Gibson et al. [ GDL93b GGM+01 ] and is in 
virtually perfect agreement with (^7|). 

For detachment into a p-wave, the general formula (pO[) reduces to the expressions: 



Jio{E) - 



M{PFf 



[2e2 Ai(e)2 -4Ai(e) Ai'(e) - 2eAi'{ef 



(88) 



M{PFf 



[2e" Ai(e)2 - Ai(e) Ai'(e) - 2e k^{ef] 



(89) 



which differ merely in a si ngle p r efactor T hese currents were first obtained via Fermi's golden rule (|1( 
by Slonim and Dalidchik |SD76, MFSFOO], while integral representations were found by Du and Delos 
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Figure 2: Some photocurrent density profiles in p-wave detachment. Clockwise, from left upper image: 
Photoelectron distribution for parallel or 7r-polarisation of the laser beam (|3|), expected distribution for 
perpendicular or cr-polarization profile for circular polarisation (|5|), and linear polarization, tilted under 
45° from the direction of field (|8q). The area displayed in each image is 1.2 mm x 1.2 mm. Parameters: 
Source-screen distance z = 0.514 m, initial electronic energy E = 60.8 /ieV, electric field strength E — 
116 V/m. 
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Figure 3: Total current (in arbitrary units) for p-wave photodetachment as a function of the electron excess 
energy E, given in ^eV, in an external electric field E = 116 V/m. Bold line: Current Jio(-E) (^8|), observed 
in TT-polarisation of the laser beam; light line: Ji.±i (|89|), valid for cr-polarization. For both the circular 
and tilted polarizations (Js^, the averaged current is predicted (dotted hne). 



[DD88, DD89] and Gibson et al. [3DL93a|. Within the experimental error margins, both expressions again 
reproduce the experimental p-wave detachment spectrum [GDL93al. 

Due to the cylindrical symmetry of the potential, the total photocurrent J{E) for arbitrary orientation of 
the polarization vector e can be expressed as a linear combination of ( |8^ ) and (|9|) (see Section 2.2). In par- 
ticular, the current Jiq{E) is observed in 7r-polarization (e || E), while Ji,±i(i?) prevails in cr-polarization 
(e _L E). In both the tilted and circular polarizations of the laser beam discussed above (p2|), the multipole 
currents contribute equal weight: Jcirc{E) = Jti\t{E) = [Jio(£') + Ji,±i(-E)]/2. For illustration, these 
currents are plotted in Figure H. While generally following their low-field limit given by the Wigner law for 
j)-wave emission J{E) cx ij"^ (p^, a strong "staircase" modulation is observed in 7r-polarization that is 
absent for the cr-case e _L E. We predicted this behaviour in Section 3.3 (^3|). 



5 Atom laser with rotating source 



In [ KBK02 ] we discuss in detail a model for an atom laser supplied by an ideal EEC. In this section we 
will summarize some of the results before we use the presented formalism to extend the theory to higher 
multipoles. For the atom laser we consider a two state model, where one magnetically trapped state (the EEC, 
acting as a quantum source) is weakly coupled to a non-trapped state (the atom laser) by a monochromatic 
radio-frequency field. (Replenishing the source condensate to maintain a truly stationary flux in the atom 



laser poses difficult technical problems. See [CSL^02].) We are interested in the propagation of coherently 



released atoms from the condensate in a homogeneous gravitational field. Explicit closed-form expressions 
for the wave function of atoms released from an isotropic Gaussian condensate containing N particles: 

cr(r) = ^/N Kla-^/^TT-^/'^ eyiY>{-r^ l{2a^)) , (90) 
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where hfl denotes the interaction strength, are given in [ KBK02 ] . The resulting beam wave function ^p{r) = 
/ d^r'G(r, r'; E) o'(r') is closely related to the wave function of the point sources discussed in Section |[ 
Similar to extended free-particle sources (38), in the far-field sector a ballistic Gaussian source may be re- 
placed by a virtual point source at a position displaced from its center along the field F. With respect to 
the virtual point source, we introduce a set of shifted variables depending on the scaled condensate width 
a — f3Fa (where /? is defined in (Ssl)), and F ~ mge^ now represents the gravitational force on a BEC atom: 



C = C + 2a^ 



(91) 



Using these variables, the wave function of the atomic beam i/)(r) (^ and the total current Jtot{E) (^ carried 
by it are convenient ly expres sed in terms of the special functions Qfc(/5, (; e) ( pl| ) and Qii.(e) ( ^9| ) discussed 
in Appendix M (cf. [KBKOl): 



(92) 



We note that the virtual source strength A(e) is strongly energy- and size-dependent: 



(93) 



(94) 



In the following we want to use the formalism of Section 2.2 to obtain closed solutions for higher states of 
angular momentum. These states naturally arise in a rotating condensate. For simplicity, we assume that the 
rotating BEC quantum fluid is in its thermodynamical ground state. This state is known to exhibit a set of 
vortices (at least, one) symmetrically arranged in an extended lattice structure [ ASRVKOl]. 



5.1 Ideal atom laser from a single vortex 



As a first example we describe the atom laser beam arising in the presence of a single vortex with fixed 
direction in a non-interacting boson gas. The wave function of the rotating condensate is then given by a first 
excited radial harmonic oscillator state [FSOl ]. (In practice, the vortex line is not stationary, but may precess 
slowly in time [HACCOl].) For simplicity, we will start with a vortex line aligned to the field direction Cz in 
an isotropic trap. In analogy to p-wave photodetachment (Section 0), the source wave function 



0-11 W 



Nhn, 



'5/2^-3/4 



^-^/^(.T + iy)exp(-r7(2a")) 



(95) 



is the oscillator eigenstate with angular momentum / — 1 and m — 1. The condensate density | cti i (r ) | drops 
to zero along the vortex line. In analogy to the multipole source formalism (Section 2.2), this wave function 
may be obtained by applying the (1,1) spherical tensor gradient K 1 1 ( V ) (|l3|) to the vortex-free BEC ground 
state (0): 



an(r) = -a\/yKn(V)a(r). 



In general, we define Gaussian multipole sources aim (r) analogous to ( 15 ) via: 

<Ji„,{r) = Ni K,„,(V)a(r) , 



(96) 



(97) 
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where A^j-^ ~ 27r'^/^a^'/r(' + 3/2) is obtained from the normahzation condition J d'^r |cr;„i(r)|^ = N{hil)'^. 
In fact, the source functions thus generated turn out to be the lowest-lying oscillator eigenstates of (I, m) 
spherical symmetry (with energy E — {I + 3/2)h'^ /Ma^). (In particular, croolr) = cr(r) covers the isotropic 
Gaussian source (^.) From (|) we evaluate the wave function ipim (r) of the corresponding outcoupled state: 



d\'a{v') Ki„^{V')G{v,v';E) 



(98) 



Here, we integrated by parts to shift the spherical tensor operator to the Green function ([4-l[). The substitution 
rule (Q) then enables us to further evaluate the integral in terms of derivatives of the known atom laser wave 
function for zero angular momentum ?/'(r) (p2|): 



^i,n{r) = Ni Ki,n [d^,dy, - Eds] ip{r) 



(99) 



Thanks to the differentiation rules ( |132| ) and (|133b, in the far-field sector these derivatives can be expressed 
as a sum of the auxiliary functions Q^.{p, C; ?) ( 130 ), analogous to the case of multipole point sources (5;m(r) 
(52). For a single vortex, it suffices to consider the sources crim(r) with I — 1, |m| = 0, 1: 



■010 (p) 



AV2P{PE)^aAii) \2C q,{p, C; i) - W Q,ip, C; i) + Qo(/5, C; e) 



1±1 



(p) = TSf3iPEfaAii)(^i±iv')Q^{pX;e), 



(100) 
(101) 



(For a = 0, point sources are recovered, and ( |100[ ), ( |101| ) become proportional to the p-wave multipole 
Green functions (^3|), (54).) If the vortex is parallel to the gravitational field, the density of the atom laser 
will be zero along the axis as it is the case for the condensate density: The vortex line is preserved in the atom 
laser profile. 

In a fashion similar to Section ^j) , we calculate the overall outcoupling rate as a function of the radiation 
frequency detuning {E = hAi'). According to (|^), the total multipole current is available from Ji„i{E) = 
— 2S[(cr;„i|G|CT;„i)]/fi.. In practice, the calculation is best carried out in momentum space, where both the 
Gaussian multipole source (^ ) and the ballistic propagator ( p^ take on a particularly simple form. Similarly 
to the point source currents (5^ (to which they reduce as a ^ 0), the Gaussian multipole currents are 
expressed using the auxiliary functions Qifc(e) (59) discussed in detail in Appendix B.2 Within the / = 1 
triplet, they explicitly read: 



Jio(e) - ^/3(/3F)3a2A(g)2 



+ Sa* Qii(e) - QiQ(g) + i Qi_^(e~) 



32 



(102) 
(103) 



Thanks to the preserved rotational symmetry of the system, all total current matrix elements Jim.i'm' (e) with 
m m' vanish, as indicated in Section 2.2. 



If £ is large, i. e., for extended condensates with a ^ 1, it is sufficient to use the asymptotic exponential 
form of the functions Qifc(e) (143). Further expanding the currents around their maximum near e = 0, we 
obtain their large-source approximations: 



a 



(104) 
(105) 
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As expected from our earlier results for a simple Gaussian source [ KBK02 1 , these currents can be interpreted 
as the integrated condensate density along a slice through the BEC at a height z fixed by the "resonance 
condition" E + Fz — (Franck-Condon principle): 



5{E + Fz) 



(106) 



(A semiclassical derivation of the "slicing approximation" is presented in Ref. |JS02|.) We note that (106) 
evidently fulfils the sum rule for the total outcoupling rate. 

For illustration, we consider two orientations of the vortex with respect to the gravitational force F. A 
vortex parallel to the field is simply represented by the Gaussian condensate wave function crii(r) (psf), and 
the ensuing laser beam characteristics are expressed in (101 ) and (103). We also examine the case of a vortex 
along the x-axis, i. e., perpendicular to F. The corresponding BEC wave function crix(r) is connected to 
the parallel vortex model by a rotation cxp(— i7rZj,/2). Application of the I ~ 1 rotation matrix for angular 



momentum eigenstates [ |3dm57 1 yields the following source term: 



1 r 



(Tii(r) + \/2(Tio(r) + cri,-i(r) 



with the associated total current: 



(107) 



(108) 



Figure ^ depicts atom laser density profiles generated by an ideal ^^Rb BEC of width a = 2 pirn at a distance 
z = 1 mm in the center of the resonance { v — 0) as well as for positive and negative detuning {v = ±4 kHz). 
For this choice of parameters, a « 3.33, so the effective energy e = ia'^ (^T|) of the assigned virtual point 



source indicates tunneling emission. According to the results presented in Section |3.4t , we expect the density 
profiles to be of overall Gaussian shape, with a mean width _D(z) = 2z/k (^s]), where z = (C + 2a'*)//3F and 
k — ApFa^ denote the distance from the virtual source and the virtual evanescent wave number, respectively. 
(Alternatively, this result follows f rom the sp reading of a minimal uncertainty wave packet of width a during 
its time of flight T = y^2Mz/F | |BBG+9^ ].) This Gaussian envelope is modulated by a factor /(^, v) that 
depends on the relative orientation of the vortex and the gravitational force. Indeed, a fairly cumbersome 
calculation yields for the asymptotic shape of the density profiles generated by the sources ( |9^ and ( 107 ) in 
the far-field sector, valid for a ^ 1: 



■ exp 



(109) 



where the modulation factors /n v) and fi± v) for paraUel and perpendicular orientation read, respec- 
tively: 



2\/2c 



(110) 



(We used the dimensionless coordinates introduced in (p4|).) Clearly, v) effects the propagation of the 

vortex in the parallel case. However, the dependence of fi±{£,, v) on the source distance ( and the detuning 
1/ = —e/2h(3 renders the atom laser profiles generated in perpendicular orientation more intriguing (see 
Figure ^: The detuning-dependent, isotropic contribution competes with a shifted parabolic term that 
grows linearly with the detector distance (. For ( <C 2a^, detuning blurs the simple vortex image present at 
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center resonance, while for C, ^ 2a'' the latter term in (IIC) dominates, causing the appearance of a node 
line in the profile whose relative position shifts linearly with the detuning i'. The transition between these 
markedly different regimes of the atom laser occurs at considerable distance from the BEC. In our example, 
we find ( — 2oL^ for a separation z « 150 /im, and the figure depicts the far-field behaviour We note, 
however, that this characteristic distance grows with the fourth power of the source size a and quickly reaches 
macroscopic dimensions: For an atom laser supplied by a BEC of width a = 10 /im, we find z « 10 cm! 

We now turn to th e frequency dependence of the total outcoupling rate. Since a ^ 1, use of the asymp- 
totic description ( 106 ) is in order In Figure |] we show the resulting current characteristics for both condensate 
orientations. While in parallel orientation the current distribution is simply Gaussian, it features a dip in the 
total current at A t^ — for a vortex line perpendicular to the gravitational field. This behaviour is easily 
understood from (106): Due to the presence of the vortex line in the slicing plane, at z = the condensate 
density adopts a minimum. We finally note here that the slici ng ap proxi mati on fails for small condensates 
(with a ^ X) and it becomes necessary to use the exact results (102) and (103). The transition between both 
regimes is studied in detail in [KBK02). 



5.2 Vortex lattices 



Using the tools developed in the previous section, we proceed to give an example of a multipole source where 
larger values of the angular momentum are present. Rotating Bose-Einstein condensates show superfluid 



behaviour and respond to an externally induced rotation by the formation of a vortex lattice [VIAH+99 



MCWDOOj , [A^SRVK01[ |ECHC02|. (A review on vortices in Bose-Einstein condensates is given in [ [FS01| ].) 
No attempt at a description of the formation and parameters of this lattice will be made. Rather, we present 
a theoretical model for a stationary atom laser supplied by an already formed vortex lattice, with vortex lines 
oriented parallel to the gravitational field F. The wave function of the lattice state (the laser source) is most 
conveniently modeled in the rotating reference frame (rf), where it becomes time-independent; we denote it 
by CTiatt.rf (r). However, the laser is observed in the laboratory frame and therefore we first have to discuss 
the transformation between both frames. 



5.2.1 Transformation to the laboratory frame 

In the laboratory frame (If), the stationary source function for the rotating frame uiatt.rf (r) becomes explicitly 
time dependent. The transformation between both frames of reference involves a uniform rotation around the 
z-axis. This rotation is generated by the unitary operator exp (— ii^i^roti/fi-), where is the z component 
of the angular momentum operator and fiiot is the frequency of rotation. The full time-dependent source term 
in the lab frame consequently reads: 

CTiatt,if(r,t) = exp(-iM/;i)exp(-iL2f2i.oti/^)o'iatt,rf(r) • (HI) 

The laser wave function resulting from a time-dependent source in the presence of the gravitational field is 
generated by the ballistic propagator (|4^: 

V'iatt,if(r,i) =-\\ d'r'y" di'if(r,i|r',i')^iatt,if(r',i') . (112) 

Introducing the time of flight T — t ~ t', and noting that the generator of the rotation commutes with the 
Hamiltonian of the accelerated motion i/grav = p^/2A/ — Mgz, this may be rewritten: 

V'iatt,if(r,i) = _ie-i(^+^^^^-)*/''^°"dre^(^+^^"-)^/'' J d^r' /^(r, r|r', 0)aiatt,.f (r') . (113) 
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Figure 4: Atom laser density profiles for a rotating ^^Rb BEC source sustaining one vortex. Left column: 
Parallel orientation of vortex li ne an d force F, with source crii(r) (p^. Right column: Vortex perpendic- 
ular to F, as given by <Ti±{r) (107). The detuning frequencies are —4 kHz (top row), kHz (center 
row), and +4 kHz (bottom row), respectively. The brightest spots of the distribution pertain to a density of 
2.5 atoms//im'^. Displayed area: 30 /im x 30 /im, distance from source: z = 1 mm; source parameters: 
a = 2 ^J.mM = 2^T X 100 Hz, N = 10^ atoms. 
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Figure 5: Total current generated per atom in a ^^Rb BEC with one vortex as a function of the detuning 
frequency /S.v. SoHd line: The source (7ii(r) (95) with the vortex line parallel to the gravitational field 
F. Dashed line: Vortex line perpendicular to F, represented by (Ti^(r) (107). Parameters: a ~ 2 fim, 
n^2Trx 100 Hz. 



Thus, the outcoupled laser wave function V'latt.if (r, t) is again stationary in the rotating frame: Like the BEC, 
the atomic beam profile rotates uniformly with frequency iliot- 

In the next step, we decompose criatt.if (r) into a superposition of angular momentum eigenstates: 



CTlatt,rf(r) = ^Cr„i,rf(r) , 



(114) 



where cTm.iflr) = mh(T,n.ii{j^)- In the laboratory frame, the rotating source function appears split into 
magnetic sublevels shifted in energy: 



(115) 



Employing the propagator representation of the ballistic Green function (^, we obtain for the beam wave 
function: 

^lattjf (r, t) = ^ ^~i{E+mm^^,)t/n f ^3^/ q^^^ mh^rot) <J^M^') . (1 16) 

The rotating source function thus yields to a description in terms of stationary sources, and the results of 
Section ID apply. 



5.2.2 The source function in tlie rotating frame 

The vortex state of the BEC is commonly described as a superposition of angular momentum eigenstates 
of the harmonic oscillator |BR99 ^oOl ]. The number of vortices and their positions are available from 
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imnimizing the energy functional in the rotating frame. For a parallel arrangement of vortices and field F, we 
may model the vortex state as a product of a two-dimensional "lattice function" (J2d{x: y) detailing the vortex 
positions (xk, yk) with a Gaussian envelope enforced by the harmonic trap potential. Introducing complex 
coefficients Vk = Xk + iyk, the lattice function is obtained as a product involving all vortex positions that 
alternatively may be expressed as a polynomial in {x + iy), 

n n 
k=l k=0 

where the coefficients w^"^ are connected to the vortex positions by the recursion relation w^"^ = 
"u^t^-i^ ~ Vk+iW^k^^K where w'"^^ = 1. (Usually, these lattices show rotational symmetry which enforces 

in) 

selection rules on the , leaving only few nonvanishing coefficients.) The complete three-dimensional 
source function in the rotating frame then reads: 

o-iatt,rf(r) = A^«exp ( — — j a2D{x,y) ■ (118) 

The constant Nn is determined by the normalization condition J d'^r |o'(r)iatt,rf P — N{hn)'^: 



, („9, 

Equation (117) then yields the decomposition of criatt,if (r) into eigenstates crm,rf (r) of ( 1 14): 

/ 2 2 2 \ 

a™,rf(r) = N^wi^iHx + iy)" exp [-^^^ - ■ (120) 

Thus, the n + 1 source components are all eigenstates of the harmonic trap potential, and the highest quantum 
number m equals the number of vortices present in the BEC. 

Here, we examine the special case of an isotropic trap (gx = = a), for which the theory outlined in 
the preceding section will deliver the outcoupling rate as well as the beam profile in analytic form. To this 



end, we note that the source components (im.if (r) (120) then simultaneously present eigenstates of and 



the total angular momentum L with quantum number I — m. Hence, the source is entirely made up from 
circular Gaussian multipole states (Tmm(r) = Kmm(V)cr(r) (pO|), (|97|): 



According to ( 1 16 ), the rotating beam is thus produced by a weighed superposition of the stationary sources 
<7mm(r), where the effective energy E„i = E + mhfliot of the various components depends on their quantum 
number m. As explained in the introduction to Section ^ outside the source region each Gaussian multipole 
source crim(r) ( p7[ ) may be mapped onto a corresponding displaced virtual point source of adjusted strength 
A(e). This allows to calculate the wave function i/'mm(r) generated by crmm(r) along the lines presented 



Gmm{r,o;E) (|52|) 



in Section [3.2| , and the final result closely resembles the corresponding balhstic multipole Green function 
): 

^mAr) = -Af3iPFfAii^) [2a + iu)]" Q„,+i(p, C; Im) , (122) 
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where e,„ = —2(3 Em + 4a'' 
from a calculation in the spirit of Section BJ 



)T|). Similarly, the total cuiTent JmmiEm) generated by CTmm(i") is available 



Jrmn (-^m ) 



n 



i(em) ■ 



(123) 



(For m = 0, 1, these expressions reduce to the results (|92[)-(P3[) and (101), (103) presented above.) Substi- 
tuting (121) and ( 1 22 ) in ( 1 1 6 ), the wave function of the rotating atom laser beam ultimately reads; 



■01att,lf (r, ^) = 51 ^ 
m=0 



-i(E+m?in,ot)t/fi. 



I (") n 

to! Wm a 



(124) 



Since cylindrical symmetry enforces that all off-diagonal elements of the total current matrix Jim,i'm'{E) 
(23) vanish (see Section 2.2), the (stationary) outcoupling rate J\a.tt{E) reduces to a properly weighed sum 
of the ballistic multipole currents Jmm{E,n) (123): 



-^latt(-E) 



1 

t V 



(125) 



We illustrate these results using a model condensate featuring a symmetrical triangular lattice of 37 vor- 
tices 10 fim apart, embedded into a Gaussian source of width a — 5 fj,m. While the frequency dependence 
of the outcoupling rate merely shows the familiar Gaussian character (compare Figures |] and plots of the 
resulting atom laser profile exhibit rich detail (Figure Due to the rotation of the source, the vortex pattern, 
which is fully transferred from the BEC into the laser beam, forms an intertwined braid-like structure along 
the z-axis. It modulates the lateral beam profile which now strongly depends on the detuning frequ ency 
Ai^. The outcoupling rate varies between the different angular momentum components (imaf (r) (121) that 
make up the source, as explained in Section 3^ ({75|). A negative shift in the frequency suppresses states with 
high |to|, leading to an approximate Gaussian shape of the particle distribution, whereas positive detuning 
(Av > 0) emphasizes these contributions. The centrifugal barrier effective for them then produces a ring-like 
"crown" emission pattern. 



6 Conclusion 

In our study, we established a systematic approach to scattering processes that involve non-isotropic emis- 
sion of quantum particles. Starting from the stationary Schrodinger equation (^ that incorporates a source 
term cr{r) responsible for a steady particle flow, we proceeded by analogy with potential theory and intro- 
duced point-like "multipole sources" (5;m(r) as limiting cases of sources with {I, to) orbital symmetry. These 
sources, and the scattering waves and currents generated by them, are available from the conventional Dirac 
^-singularity and its assigned Green function G(r, r'; E) through the action of a differentiation operator 



of suitable spherical symmetry, the spherical tensor gradient Kim{^') (Section 2.2). The scattering waves 



emitted by multipole sources locally show pure (Z, to) angular symmetry, and in the absence of an external 



potential, they reduce to the spherical partial waves familiar from conventional scattering theory (Section 2.4 ). 

Notably, the local orbital characteristics are preserved also in an external potential that breaks rotational 
symmetry. The multipole wave then describes the propagation of particles initially emitted in the (Z, to) 
eigenstate of angular momentum, and thus generalize the notion of a partial wave. We performed a detailed 
study of the linear potential environment U{y) — — r • F, i. e., scattering in the presence of a uniform 
force field F. This problem allows for an analytical solution, and closed-form expressions for the ballistic 
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Detuning Frequency Az/ = E/h [kHz] 



Figure 6: Total current per atom in a *^Rb BEC with 37 vortices as a function of detuning frequency Ai/. 
The vortex structure of the BEC is not visible in the integrated current. Vortex separation in the BEC: 10 /xm, 
rotation frequency 57iot = 27r x 250 Hz, outcoupling strength f2 = 27r x 100 Hz, size parameter a = 5 ^m. 




Figure 7: Beam profile of an atom laser from a BEC with 37 vortices at different energies and heights. 
From the left figure to the right one we vary the detuning frequency (—10, 5, 20) kHz. These frequencies 
are also marked by dotted lines in Figure |[ The density in the beam is plotted at three different distances 
(177, 239, 300) /^m from the center of the BEC. The vortex cores are indicated by small tubes at all distances. 
BEC parameters same as Figure |^. 
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multipole waves Gim{j^,o; E) and cuiTents Jim,i'm'{E) ( |60| ) are assembled in Section |[ Uniformly 
accelerated scattering waves display a characteristic set of features, including a prominent fringe structure and 
a modulation of the cross section, that in a semiclassical explanation are attributed to two-path interference 



in the force field (Section 3.4) 



The theory of ballistic multipole waves directly applies to near-threshold photodetachment processes 
in an electric field environment, a topic that attracted considerable interest during the past fifteen years. 
The source model immediately yields analytical expressions for the photocurrent spectrum and the spatial 
electron distribution that are in excellent agreement with the available experimental data. Our attention was 
mainly dedicated to p-wave photodetachment where we discussed the influence of the laser polarization on 
the photoelectron current profile (Section 

Somewhat akin to the situation in potential theory, in the ballistic environment extended sources of Gaus- 
sian shape may be replaced by "virtual" point sources that are displaced from the center of the actual dis- 
tribution. This mapping remains feasible for "Gaussian multipole sources" cr;m(r) (p7|), harmonic oscillator 
eigenstates of (/, m) orbital symmetry that are generated by the same differentiation formalism as the multi- 
pole point sources 6im{r). Depending on their size, these Gaussian multipoles may act as effective ballistic 



tunneling sources that are characterized by very different properties (Section p.'% . In practice, Gaussian wave 
functions are shared by the atoms in an ideal Bose-Einstein condensate trapped in a harmonic potential. A 
ballistic tunneling source is then realized by continuous outcoupling of atoms under the influence of gravity, 
leading to the formation of an "atom laser" beam. While a BEC in its ground state simply leads to isotropic 
emission, vortices embedded in the condensate will create higher angular momentum modes in the laser pro- 



file. In particular, a BEC bearing a single vortex provides a p-wave ballistic tunneling source (Section 5.1). 
The particle distribution generated by it depends strongly on the relative orientation of vortex and gravita- 
tional force: In parallel alignment, the vortex is simply preserved in the profile, while the perpendicular setup 
is characterized by an unusual detuning-dependent node structure in the atom distribution. Finally, we also 
discussed the properties of an atom laser outcoupled from a rapidly rotating BEC sustaining a vortex lattice in 



Section 5.2. The rotating beam wave function thus created is a coherent superposition of ballistic multipole 
waves with very high angular momentum, and the resulting laser profile starkly depends on the detuning from 
resonance, while the total outcoupling rate shows little variation. 
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A Translation theorem for harmonic polynomials 

In this appendix, we present a formula that allows to expand a harmonic polynomial Kim (r + a) with shifted 
center into a spherical power series with respect to the origin, i. e., in the variable r. Since the position 
variables r and a in the argument are interchangeable, the same series will also furnish the expansion of 
Kim(r + a) around a. Hence, we expect the general form for this series: 

; A 

Kimir + a) = ^ Kam(i-) Ki_;,,,„_^(a) , (126) 

A=0 p=-A 

where Cj™ = C^i^x m~u ^^^^t hold. (Note that Kj^ (r) is a homogeneous polynomial of order /, so the orders 



in the right-hand side products of harmonic polynomials in (126) must add up to /. Similarly, the sum of their 
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magnetic quantum numbers mu st be to.) The coefficients Cj^™ in (126) have been established in closed form 



inRefs. [Cao78, CD95, Bra99]: 



Clrn 



47r(2Z + 1) 



I — TO 



(2A + 1){21 -2X+1)\X + ^IJ \\- 



(127) 



In particular, Cqq ~ V4tt. (For |to — /i| > / — A, the coefficient vanishes.) 

The general series (126) simplifies if the shift in the argument of K;m(r + a) takes place along the axis 
of quantization, i. e., a = acz- Then, rotational symmetry around the z-axis is preserved, and the quantum 



number to is not affected by the translation. Hence, only terms with /i = m survive in (126). Inserting the 
explicit value KAo(ae^) = ^/i2X + l)/47ra^ [ |Mes64| ], we thus obtain [ ^ob31[ , [VIF53p : 



K;„(r + aiz 



"■'■K,m(r) , 



where the translation coefficient Tjim is given by (127): 



21 -2j + l 



47r ]j 2j 

These developments prove useful throughout Section ||. 



2/ + 1 (l + m\ fl-m 
l\j + mj\j — m 



(128) 



(129) 



B Some integrals involving Airy functions 



B.l The functions Qfc(p, C; e) 

In the course of our investigation into ballistic multipole matter waves, integral expressions of the following 
type: 

dr 



Qfc(P>C;e) 



exp < 1 



p2 r3' 



(130) 



are frequently encountered. For integer indices k, this set of integrals permits explicit evaluation in terms of 



products of Airy functions. Introducing the Airy Hankel function Ci(M) = Bi(u) + i Ai(M) [ AS65|, the basic 
member of this class reads: 

Qo(p, C; e) = Ai(e - C + p) Ci(e - C - p). (131) 

This result is easily verified by observing that ( |130| ) in this case reduces to the Laplace transform of the 
ballistic propagator in one dimension (using dimensionless units) and therefore must equal the well-known 
Green functi on of a freely falling particle in a single spatial dimension [ Rob96 1 . Indeed, apart from scaling 
the integrals (130) for positive integer indices represent the stationary ballistic Green functions in the spaces 
of odd dimension D — 2fc + 1. In particular, uniformly accelerated waves in physical space (D ~ 3) are 
represented by the function (p, (; e) [ |DS76| , [BBG+98| , |KBK02[ ]. 

From the definition ( |130| ), two recurrence formulae for increasing and decreasing value of the index k are 
immediately available: 



Q/c+i(p>C;e) 



2_d_ 
Yp 'dp 



Qfc(p,C;e) 



l_ d_ 
Yp 'dp 



(132) 
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Q-(fc+i)(Pj ~ 9( 



Qk 



-A;(p,C;e) = Qo(p>C;e)- 



(133) 



(The latter equalities require k > 0.) From a practical point of view, the expressions thus obtained become 
rather unwieldy with growing \ k\. The following five-point recursion relation, again easily verified using the 
integral representation (130), presents a favorable alternative: 



Qfc+2(p, C; e) - (fc + Qk+M C; e) + (C - e) QJp, C; e) + J Qk-2{p, C; e) 



0. 



(134) 



(In p assing, we remark that the general solution to ( 134 ) is given by the recurrence relations ( 132 ) and 
( 133 ) once we replace the basic function Qo(/3, C; e) (131) by the more general expression Qo(/5, C', e) = 
Zi(e — C + p) 3i(e — C ^ p)> where Zi{u), 3i(") denote two arbitrary solutions of the Airy differential 
equation, i. e., linear combinations of the Airy functions AUu) and Bi{u).) 

Finally, we inquire into the asymptotic behaviour of (130) in the limit p ^ 0. Here, we are interested in 
the case of integer index fc > 1. Then, small values of r provide the bulk contribution to the integral, which 
allows us to neglect the linear and cubic terms in the exponent of (130). In this approximation, the integral 
evaluates to a Gamma function of half-integer argument [ AS65 1 : 



Qfc(p,C;e) 



r(fc-l/2) (2fc-3)!! 



(135) 



Therefore, the function Qj,(/3, C; diverges as p ^ 0. This singularity, however, affects only the real part of 

Qfc(p,C;e)- 



B.2 The functions Qi^. (e) 

Another i mpo rtant class of functions that regularly appears when calculating ballistic total currents is con- 
tained in (130) as a limiting case: 



lim lim 5{Qfe(/3,C;e)}. 



(136) 



Unlike the functions Qfc(/0, C; e) that are divergent in this limit for fc > (reflecting the multipole source 
singularity), their imaginary parts Qife(e) remain well-defined. Obviously, Qio(e) = Ai(e)2, and all other 
expressions are available from suitably modified recurrences (132), (133): 



1. A 

2i 



Ai(e + z) Ai(e - z) 



(137) 



lim { —r Ai(e - zf 



(138) 



(fc > 0). Thus, the functions Qife(e) can be extracted from the Taylor series of Ai(e + z) Ai(e — z) and 
Ai(e — z)2, respectively. For practical purposes, again the following recursion relation adapted from (134) 
proves more suitable: 



Qife+i(e) + eQife(e) 



^ife-2(e) = 0. 



(139) 



37 



For the sake of completeness, we note that unhke the functions Qj.(p, (; e) (13C), the limits ( 136 ) can also be 
evaluated for half-integer index Qife+i/2(e), which in turn allow s to c alculate ballistic multipole currents in 
spaces of even dimension. Here, direct evaluation of the integral (OQ) shows that: 



Qii/2(e) 



1 



1 



Aii(22/3e) 



(140) 



holds, where Aii(w) = Jq dz Ai{z) denotes the integral of the Airy function [AS65|. (In particular, 
Qii/2(0) = ^/6^/^^.) The other functions of half-integer index Qifc+i/2(e) are available from the differ- 
entiation formula: 



Qi, 



fc-l/2 



d_ 



k+l/2 



(141) 



as well as the recursion relation (139) 



B.3 Approximations for large values of e 



Next, we aim to establish the asymptotic behavior of the functions Qife(e) for integer k in the limit |e| oo. 
The expressions thus obtained are helpful in the discussion of the ballistic multipole currents presented in 
section |^. We start out with the observation that the integral representation of these functions in ( 130 ) and 
(|136[) may be rewritten in a complex contour integral: 



1 



du 



C 



,k+l/2 



(142) 



Here, the paths C leads along the imaginary axis, avoiding the singularity at the origin and the cut in the 
complex plane which we choose to place onto the negative real axis (see Figure ||). (Note that the integrand 
in (142) is double-valued. This complication is absent for half-integer indices.) 

For large values of |e|, a saddle point approximation for (142) is in order. This poses no problems for 
e +CX) (tunneling case). Then, the relevant stationary point of the exponent S- = l^fl is readily identified, 
and shifting the path of integration (see Figure ^ ultimately yields the asymptotic series: 

( \ 1 ^\-(fe+l) ( 4 3/2\ f- 3fc2 



1 



9fc + 5 



24e3/2 



O 



(143) 



The situation is more involved for e — > — oo (classically allowed motion). Here, the saddle points are located 
at iSi.iii = =F2i\/T^' and the direction of steepest descent cuts the imaginary axis under an angle of 7r/4. 
Hence, the integration path must be deformed to lead into the sector of asymptotically vanishing integrand 
I argwl > 57r/6. Unfortunately, due to the presence of the cut in the complex plane, the partial paths Ci, Cm 
cannot be simply connected as di[u] —oo. Rather, they must be linked by an additional path element Cu 
that loops back around the singularity located at w — 0, as indi cated i n Figure |[ The latter contribution is 



asymptotically evaluated by means of Hankel's integral formula [ ]01v74 ] that states J^,^^ e*t ^dt = 27ri/r(2;): 



du 



Cii 



,k+l/2 



-isu+«''/12 



27ri 



|fc-l/2 



r(fc + i/2)' 



(144) 



This secular part obeys a simple power law dependence reminiscent of the Wigner law (2.4, [ Wig48]), while 
the saddle points 5*1, 5*111 deliver contributions of oscillatory character: 



Cii: 



4i 



|3/2 



(145) 
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Figure 8: Evaluation of the integral (142). The figure displays the original contour C and its shifted coun- 



terparts for large values of e. For e > 0, the contour is deformed to C_ and runs through the saddle point 
S- = 1\ft- In the case e < 0, the contour is split into three parts Ci, C\\, Cm traversing the saddle points 
S\, 5*111 at ±2iY^, yet avoiding the cut in the complex u plane (dark grey line) by circling the singularity at 
M = 0. The integrand asymptotically vanishes in the shaded sectors. 



(cf. eg. ( 143)). The saddle point S\ adds the conjugate complex result, and upon gathering the terms in (144) 
and ( |145| ), one obtains the following leading asymptotic form for Qifc(e) as e ^ -co: 



|fc-l/2 



20F r(fc + l/2) 



27r V 



-(fc+i) 



sm 
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References 

[AS65] M. Abramowitz and LA. Stegun. Handbook of Mathematical Functions. Dover, New York, 
1965. 

[ASRVKOl] J.R. Abo-Shaeer, C. Raman, J.M. Vogels, and W. Ketterle. Observation of vortex lattices in 
Bose-Einstein condensates. Science, 292:476, 2001. 

[Bay78] B.F Bayman. A generalization of the spherical harmonic gradient formula. J. Math. Phys., 
19:2558, 1978. 

[BBG+98] C. Bracher, W. Becker, S.A. Gurvitz, M. Kleber, and M.S. Marinov. Three-dimensional tunnel- 
ing in quantum ballistic motion. Am. J. Phys., 66:38, 1998. 

[BBW73] T. Banks, CM. Bender, and T.T. Wu. Coupled anharmonic oscillators. I. Equal-mass case. 
Phys. Rev. D, 8:3346, 1973. 



39 



[BDD96] C. Blondel, C. Delsart, and F. Dulieu. The photodetachment microscope. Phys. Rev. Lett, 
77:3755, 1996. 

[BDDV99] C. Blondel, C. Delsart, F. Dulieu, and C. Valli. Photodetachment microscopy of O" . Eur. Phys. 
J. D, 5:207, 1999. 

[BHE99] I. Bloch, T.W. Hiinsch, and T. EssUnger. Atom laser with a cw output coupler. Phys. Rev. Lett., 
82:3008, 1999. 

[BMS+87] H.C. Bryant, A. Mohagheghi, J.E. Stewart, J.B. Donahue, C.R. Quick, R.A. Reader, V. Yuan, 
C.R. Hummer, W.W. Smith, S. Cohen, W.P. Reinhardt, and L. Overman. Observation of 
motional-field-induced ripples in the photodetachment cross section of H~. Phys. Rev. Lett, 
58:2412, 1987. 

[BR99] D.A. Butts and D.S. Rokhsar. Predicted signatures of rotating Bose-Einstein condensates. Na- 
ture, 397:327, 1999. 

[Bra99] C. Bracher. Quantum Ballistic Motion and its Applications. PhD thesis, Technische Universitat 
Miinchen, 1999. Unpublished. 

[BRK97] C. Bracher, M. Riza, and M. Kleber. Propagator theory in scanning tunneling microscopy. Phys. 
Rev. B, 56:7704, 1997. 

[BV71] V.L. Bakhrakh and S.l. Vetchinkin. Green's functions of the Schrodinger equation for the sim- 
plest systems. Theoret. Math. Phys., 6:283, 1971. [Teoret. Mat. Fiz. 6, 392-402 (1971)]. 

[Cao78] M.J. Caola. Solid harmonics and their addition theorems. J. Phys. A: Math. Gen., 1 1:L23, 1978. 

[CD95] S. Chakrabarti and D.P. Dewangan. Addition theorems for sohd harmonics and the second Bom 
amplitudes. /. Phys. B: At Mol. Opt. Phys., 28:L769, 1995. 

[CSL+02] A.P. Chikkatur, Y. Shin, A.E. Leanhardt, D. Kielpinski, E. Tsikata, T.L. Gustavson, D.E. 

Pritchard, and W. Ketterle. A continuous source of Bose-Einstein condensed atoms. Science, 
296:2193, 2002. 

[DD88] M.L. Du and J.B. Delos. Photodetachment of H" in an electric field. Phys. Rev A, 38:5609, 

1988. 

[DD89] M.L. Du and J.B. Delos. Effect of an electric field on the photodetachment of H~. Phys. Lett. 
A, 134:476, 1989. 

[DK082] Yu.N. Demkov, V.D. Kondratovich, and V.N. Ostrovskii. Interference of electrons resulting 
from the photoionization of an atom in an electric field. JETP Lett., 34:403, 1982. [Pis'ma Zh. 
Eksp. Teor. Fiz. 34, 425^27 (1981)]. 

[DMM75] V.V. Dodonov, l.A. Malkin, and V. Man'ko. The Green function of the stationary Schrodinger 
equation for a particle in a uniform magnetic field. Phys. Lett. A, 51:133, 1975. 

[DS76] F.I. Dalidchik and V.Z. Slonim. Strong exchange interaction effects in a homogeneous electric 
field. Sov Phys. JETP, 43:25, 1976. [Zh. Eksp. Teor. Fiz. 70, 47-60 (1976)]. 

[Du89] M.L. Du. Oscillation of electron flux in photodetachment of H~ in an electric field. Phys. Rev. 
A, 40:4983, 1989. 



40 



[ECHC02] P. Engels, I. Coddington, P. C. Haljan, and E. A. Cornell. Nonequilibrium effects of anisotropic 
compression applied to vortex lattices in bose-einstein condensates. Phys. Rev. Lett., 89: 100403, 
2002. 



[Eco83] E.N. Economou. Green's Functions in Quantum Physics (Solid-State Sciences 7). Springer, 
Berlin, 1983. 

[Edm57] A.R. Edmonds. Angular Momentum in Quantum Mechanics. Princeton University Press, 
Princeton, 1957. 

[Fab81] I.I. Fabrikant. Interference effects in photodetachment and photoionization of atoms in a homo- 
geneous electric field. Sov. Phys. JETP, 52:1045, 1981. [Zh. Eksp. Teor. Fiz. 79, 2070-2077 
(1980)]. 

[Fab82] I.I. Fabrikant. Scattering of electrons by atoms and photodetachment of negative ions in a 
homogeneous electric field. Sov. Phys. JETP, 56:967, 1982. [Zh. Eksp. Teor. Fiz. 83, 1675- 
1684 (1982)]. 

[Fab91] I.l. Fabrikant. Near-threshold photodetachment of H~ in parallel and crossed electric and mag- 
netic fields. Phys. Rev. A, 43:258, 1991. 

[Fab94] I.I. Fabrikant. Photodetachment from the S~ and Cl~ ions in a static electric field and low- 
energy electron scattering by S and CI atoms. J. Phys. B: At. Mol. Opt. Phys., 27:4545, 1994. 

[FH65] R.P. Feynman and A.R. Hibbs. Quantum Mechanics and Path Integrals. McGraw-Hill, New 
York, 1965. 

[FSOl] A.L. Fetter and A.A. Svidzinsky. Rotating vortex lattice in a Bose-Einstein condensate trapped 
in combined quadratic and quartic radial potentials. / Phys.: Condens. Matter, 13:R135, 2001. 

[Gav92] M. Gavrila, editor. Atoms in Intense Fields. Advances in Atomic, Molecular, and Optical 
Physics, Supplement. Academic Press, London, 1992. 

[GDL93a] N.D. Gibson, B.J. Davies, and D.J. Larson, p-wave photodetachment in a static electric field. 

Phys. Rev. A, 48:310, 1993. 

[GDL93b] N.D. Gibson, B.J. Davies, and D.J. Larson, s-wave photodetachment in a static electric field. 
Phys. Rev A, 47:1946, 1993. 

[GGM+01] N.D. Gibson, M.D. Gasda, K.A. Moore, D.A. Zawistowski, and C.W. Walter, s-wave photode- 
tachment from S~ ions in a static electric field. Phys. Rev. A, 64:061403-1, 2001. 

[GKK91] B. Gottheb, M. Kleber, and J. Krause. Turmehng from a 3-diniensional quantum well in an 
electric field: an analytic solution. Z. Phys. A - Hadrons and Nuclei, 339:201, 1991. 

[Gol97] PA. Golovinskii. Interference in the photodecomposition of negative atomic hydrogen ions in 
an electric field. Sov. Phys. JETP, 85:857, 1997. [Zh. Eksp. Teor. Fiz. 112, 1574-1583 (1997)]. 

[Gou72] G. Gountaroulis. Green-function of the free electron in a uniform magnetic field. Phys. Lett. A, 
40:132, 1972. 



41 



[HACCOl] 

[HFCK90] 

[HM93] 
[HoOl] 

[Hob31] 

[Hor63] 
[Hos63] 

[HP64] 

[HS52] 

[HS67] 

[JS02] 

[KBKOl] 
[KBK02] 
[KO90a] 

[KO90b] 

[KP37] 
[LKB80] 



P.C. Haljan, B.P. Anderson, I. Coddington, and E.A. Cornell. Use of surface-wave spec- 
troscopy to characterize tilt modes of a vortex in a Bose-Einstein condensate. Phys. Rev. Lett, 
86(14):2922, April 2001. 

Z.H. Huang, T.E. Feuchtwang, P.H. Cutler, and E. Kazes. Wentzel-Kramers-Brillouin method 
in multidimensional tunneling. Phys. Rev. A, 41:32, 1990. 

S. Horch and R. Morin. Field emission from atomic size sources. J. Appl. Phys., 74:3652, 1993. 

T.-L. Ho. Bose-Einstein condensates with large number of vortices. Phys. Rev. Lett, 
87(6):060403-l, August 2001. 

E.W. Hobson. The Theory of Spherical and Ellipsoidal Harmonics. Cambridge University 
Press, Cambridge, 1931. 

L. Hormander. Linear Partial Differential Operators. Springer, Berlin, 1963. 

L. Hostler. Coulomb Green's functions and the Furry approximation. J. Math. Phys., 5:591, 
1963. 

L. Hostler and R.H. Pratt. Coulomb Green's function in closed form. Phys. Rev. Lett., 10:469, 
1964. 

I. Halperin and L. Schwartz. Introduction to the Theory of Distributions . University of Toronto 
Press, Toronto, 1952. 

H. M. Van Horn and E.E. Salpeter. WKB approximation in three dimensions. Phys. Rev., 
157:751, 1967. 

Y. Japha and B. Segev. Semiclassical theory of field-induced thermal transition rate with appli- 
cation to output coupling of a Bose-Einstein gas at finite temperature. Phys. Rev. A, 65:06341 1- 

I, 2002. 

T. Kramer, C. Bracher, and M. Kleber. Four-path interference and uncertainty principle in 
photodetachment microscopy. Europhys. Lett, 56:471, 2001. 

T. Kramer, C. Bracher, and M. Kleber. Matter waves from qua ntum sources in a f orce field. J. 
Phys. A: Math. Gen., 2002. (accepted for publication, e-print: quant-ph/01 10155 ). 



V.D. Kondratovich and V.N. Ostrovskii. Resonance and interference phenomena in the pho- 
toionization of a hydrogen atom in a uniform electric field: III. Comparison with recent experi- 
mental and theoretical results. J. Phys. B: At Mol. Opt Phys., 23:21, 1990. 

V.D. Kondratovich and V.N. Ostrovskii. Resonance and interference phenomena in the pho- 
toionization of a hydrogen atom in a uniform electric field: IV. Differential cross sections. J. 
Phys. B: At Mol. Opt Phys., 23:3785, 1990. 

PL. Kapur and R. Peierls. Penetration into potential barriers in several dimensions. Proc. Roy. 
Soc. A, 163:606, 1937. 

E. Luc-Koenig and A. Bachelier Systematic theoretical study of the Stark spectrum of atomic 
hydrogen II: density of oscillator strengths. Comparison with experimental absorption spectra 
in solid-state and atomic physics. J. Phys. B: Atom. Molec. Phys., 13:1769, 1980. 



42 



[LKKB97] A. Lohr, M. Kleber, R. Kopold, and W. Becker. Above-threshold ionization in the tunneling 
regime. Phys. Rev. A, 55:R4003, 1997. 

[LLF90] Y.L. Li, C.H. Liu, and S.J. Franke. Three-dimensional Green's function for wave propagation in 
a linearly inhomogeneous medium — ^the exact analytic solution. J. Acoust. Soc. Am., 87:2285, 
1990. 

[MAH+99] M.R. Matthews, B.R Anderson, RC. Haljan, D.S. Hall, C.E. Wieman, and E.A. Cornell. Vor- 
tices in a Bose-Einstein condensate. Phys. Rev. Lett., 83:2498, 1999. 

[MAK+97] M.-O. Mewes, M.R. Andrews, D.M. Kum, D.S. Durfee, C.G. Townsend, and W. Ketterle. Out- 
put coupler for Bose-Einstein condensed atoms. Phys. Rev. Lett., 78:582, 1997. 

[MCWDOO] K.W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard. Vortex formation in a stirred Bose- 
Einstein condensate. Phys. Rev. Lett., 84:806, 2000. 

[Mes64] A. Messiah. Quantum Mechanics, volume 1. North-Holland, Amsterdam, 1964. 

[MF53] RM. Morse and H. Feshbach. Methods of Theoretical Physics, volume 2. McGraw-Hill, New 
York, 1953. 

[MFSFOO] N.L. Manakov, M.V Frolov, A.F. Starace, and I.I. Fabrikant. Interaction of laser radiation with 
a negative ion in the presence of a strong static field. J. Phys. B: At. Mol. Opt. Phys., 33:R141, 
2000. 

[MJ59] G. MoUenstedt and C. Jonsson. Elektronen-Mehrfachinterferenzen an regehnaKig hergestellten 
Feinspalten. Z. Phys., 155:472, 1959. 

[Mul66] C. Miiller. Spherical Harmonics (Lecture Notes in Mathematics Vol. 17). Springer, BerUn, 
1966. 

[NOV+02] C. Nicole, H.L. Offerhaus, M.J.J. Vrakking, F. Lepine, and C. Bordas. Photoionization mi- 
croscopy Phys. Rev Lett., 88:133001-1,2002. 

[01v74] F.W.J. Olver. Asymptotics and Special Functions. Academic Press, New York, 1974. 

[Rob96] R.W. Robinett. Quantum mechanical time-development operator for the uniformly accelerated 
particle. Am. J. Phys., 64:803, 1996. 

[Row78] E.G.P. Rowe. Spherical delta functions and multipole expansions. J. Math. Phys., 19:1962, 
1978. 

[Sch73] J. Schwinger. Particles, Sources, and Fields, volume 2. Addison- Wesley, 1973. 

[Sch81] L.S. Schulman. Techniques and Applications of Path Integration. Wiley, New York, 1981. 

[SD76] VZ. Slonim and F.I. Dalidchik. Impurity electroabsorption and photodestruction of negative 
ions in an electiic field. Sov Phys. JETP, 44:1081, 1976. [Zh. Eksp. Teor. Fiz. 71, 2057-2067 

(1976)]. 

[TH83] J. Tersoff and D.R. Hamann. Theory and application for the scanning tunneUng microscope. 
Phys. Rev Lett., 50:1998, 1983. 



43 



[Wig48] E.P. Wigner. On the behavior of cross sections near thresholds. Phys. Rev., 73:1002, 1948. 

[WRG88] H.Y. Wong, A.R.P. Rau, and C.H. Greene. Negative-ion photodetachment in an electric field. 
Phys. Rev. A, 37:2393, 1988. 

[WS83] E.J. Weniger and E.O. Steinbom. New representations for the spherical tensor gradient and the 
spherical delta function. J. Math. Phys., 24:2553, 1983. 



44 



